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Abstract. This paper covers several techniques of intercompar- 
ison of Cosmic Microwave Background (CMB) anisotropy exper- 
iments and models of structure formation. It presents the con- 
straints on several cosmological parameters using current CMB 
observations. 



1. Introduction 



In 1964 Penzias and Wilson 1 143 ] discovered the Cosmic Microwave Background 
(CMB) radiation as an excess instrument noise with temperature ~ 3K, coming 
from all directions in the sky. Dicke and Peebles et al. (1965) interpreted this 
as the relic radiation from the Big Bang. The discovery of this relic radiation in 



conjunction with Hubble's observation of the cosmic expansion in 1927 [ 100 1 , the 
synthesis of the light elements, the evolution of radio source counts etc, provided 
strong support to the hot Big Bang model. 

Recently the FIRAS instrument on th e CQ BE satellite has confirmed that 



this radiation has a Planck spectrum with [L33] 



T = 2.726 ± 0.010K (95% confidence). (1) 

In the standard Big Bang model thermalization of the CMB occurs at a redshift 
z ~ 10 6 at T > 10 7 K, via double Compton scattering (e + 7 — > e + 27) and 
bremsstrahlung (e + p — > e + p + 7) processes, resulting in a Planckian spectrum 
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kept along its subsequent evolution. So, the observed black body spectrum of 
the CMB radiation suggests its origin in the early universe; this is further sup- 
ported by its high degree of isotropy. The CMB photons were last scattered at 
a redshift z ~ 1100, when the expanding universe had cooled to temperatures 
T ~ 3000K, sufficient for protons to capture the free electrons producing neu- 
tral atoms. After this period of recombination the CMB photons free stream 
towards us from this 'last scattering surface', providing direct information about 
this early epoch 300,000 years after the initial singularity. Initial theoretical pre- 
dictions suggested that structure formation in the universe like galaxies, clusters 
of galaxies, etc., implied fluctuations in the CMB, due to the presence of pri- 
mordial seed structures at the time when the CMB radiation decouples from the 
matter distribution. Consequently observations of CMB radiation anisotropies 
provide a means of testing stucture formation theories. On large angular scales 
the CMB fluctuations may probe the inflationary power spectrum and therefore 
the physics o f th e first 10~ 35 s. The original detection paper by Penzias and 
Wilson 1965 [1431, set an upper limit of < 10% on any anisotropy. The initial 
predictions indicated that galaxy form ation wo uld prod uce fluctuations in the 
CMB of the order of 1 part in 10 2 Jl63|, or 10 3 Jl7|, |l74| . 

Since then work on the CMB anisotropy question has shown improvements in 
experimental sensitivity and in theoretical calculations with increasingly sophis- 
ticated estimations of the temperature fluctuations of the CMB. For universes 



containing predominantly bary oni c ma tter, a prediction of AT/T ~ 10 was ob- 



tained by several authors 
dark matter universes 
electron neutrino mass of 




, |8b1,p0l|. By 1980, the hypothesis of non-baryonic 
was enco uraged by the claims of the existence of a 
30eV (l27j . Although the neutrino as a candidate of 
dark matter failed to explain the st ruct ure formation 103, 203, 204], it still is an 
acceptable dark matter candidate [ 181 1 . At the same time inflationary scen arios 
paid particular attention to cold dark matter dominated universes [ gjj , 141, 69[ , 
where cold means that the candidate particles have small thermal velocities (the 
free streaming length of such particles is much smaller than the size of any galaxy 
or cluster of stars) . Fo l lowing im provements in the upper limits on smaller scales 



of AT/T ~ 10" 4 [193 



191 



195 1, new theoretical calculations predicted that for 



models with a cold dark matter component, the amplitude of the temperature 
fluctuations of the CMB could be reduced by an order of magnitude. 

The first detection of large an gula r scale primordial CMB anisotropy cam e 
in 1992, reported by Smoot et al. [178) using the COBE DMR experiment [ pT?| ]. 
This initial result consisted of statistical evidence for CMB anisotropy; the CMB 
maps were not sensitive enou gh t o pin down actual CMB features on the sky. 



In 1993 the MIT experiment [ 134 , |70| reported results at high frequency which 
confirmed the COBE detection. Given the frequency range and the sky coverage 
between these two experiments this confirmed that the detected signal was in- 
deed CMB anisotropy. This was followed by low frequency confirmation in 1994, 
from the Tenerife experiment , which enabled the association of the detected 
statistical anisotropy with actual CMB features on the sky. Since the COBE 
detection, several other experiments have reported detections on intermediate 
angular scales (see Section ph, although given their sky coverage and frequency 
range spanned, the possibility of foreground contamination is in some of the cases 
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still under consideration. The systematic effects affecting some of these obser- 
vations are another cause of conc ern. Recently new i ntere sting detections have 
been reported by the Saskatoon 136] and the CAT 171] experiments, among 
others, which with the other observations cover angular scales ranging from 10° 
to 0.25°. Meanwhile questions such as the process of structure formation, the 
origin of the seed perturbations, the nature and quantity of the matter content 
of the universe, the values of the density parameter Q, of the Hubble constant Ho 
, and of other cosmological parameters remain still to be answered conclusively. 

All these new exciting results and more sophisticated theoretical calculations 
motivated the work carried out by several teams in search of satisfactory answers 
to such questions, bearing in mind the limitations imposed by the status of the 
data. 



2. Essential concepts and mathematical tools 

2.1. Statistical description of the C MB field 

It is standard procedure to express the temperature fluctuations of the CMB on 
the celestial sphere using an expansion in spherical harmonics: 

oc I 

AT/T(a, 4>) = 5^X1 a ^ Y im(a,4>) (2) 
1=2 m=— l 

where (a, <f>) are the polar coordinates of a point in the spherical surface. We have 
ignored the dipole (I = 1) which is largely dominated by the Doppler effect due to 
the motion of our Galaxy with respect to the last scattering surface. Inflationary 
models predict that the initial CMB perturbations follow a Gaussian distribution, 
and as the observed CMB field conserves the original distribution we can consider 
it as a 2-dimensional Gaussian random field on the sky. In this case the multipole 
moments ai m are stochastic variables with random phase, zero mean and vari- 
ance given by Ci = (|a! m | 2 ) which defines a power-spectrum. The properties of a 
Gaussian field are fully specified by the two-point angular auto-correlation func- 
tion (ACF), the expected average of the product of the temperature fluctuations 
in two directions in the sky separated by an angle 6 : 

I AT AT \ 
C(0) = (y(ni)y(m)) (3) 

where ni ■ ni = cos 9 and the angle brackets refer to an average over the whole 
sky. Using Equation (Q) the expression of the ACF becomes: 

C ran {9) = 2J /J (aimat'm') Y im(a, 4>)Yy m , (a, 4>) (4) 

For a temperature pattern that is statistically isotropic, the rotational symmetry 
implies: 



( a im a ;'m'} — Sll'Smm'Gt- 
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(5) 



Using the relation (JE]) and the addition theorem for spherical harmonics, one can 



express the random variable C(9), as an expansion in Legendre polynomials [167]: 

c ran (e) = ^Y. Qt( x ) p i(cos e), (6) 

i>2 

where 

i 

Qi( x ) = X! \ aim \ 2 ' ( 7 ) 

m— — / 

with a x1/h-i distribution and 

(Q?(x)) = (2Z + l)C ! , (8) 

and variance 

var(Q?) = 2{Q 2 f/(2l + 1). (9) 

The relationship between the angular scale 9 and the multipole I is given by 
I ~ 2/9, with 9 in radians. The coefficients Ci constitute the angular power 
spectrum of the temperature fluctuations and, as we will see, play an important 
role in the comparison of CMB experiments with theoretical models of structure 
formation. As shown in Equation (^) the uncertainty scales as (£ + l/2)~ 1//2 which 
explains the fact that sm aller angular scale experiments have smaller cosmic 



variances (see also Section 2.4), since they probe larger multipoles I. For large 
scale experiments this is an unavoidable uncertainty and in particular it limits 
the accuracy of the normalization of the power spectrum. The expectation value 
of the temperature correlation function, i.e., the ensemble average, is given by: 

C(0) = ^-^(2J + ljCiP^cos 0). (10) 
I 

In this statistical description of the CMB fluctuations, our universe is con- 
sidered as a realization of the theoretical field, and therefore the observed ACF 
represents an estimation of the ACF of the field from a finite sample. The un- 
certainty of these estimations constitutes an intrinsic limitation imposed by the 
statistical nature of the theory and is independent of the instrumental sensitivity. 

2.2. (AT/T) rrns 

One of the observables in CMB experiments is the rms of temperature fluctu- 
ations. This quantity may then be compared with the expected temperature 
fluctuation for a given experiment and model. For a ri gorous comparison the 



theoretical uncertainty must also be considered (see 15£ ] ) . 

The expected temperature fluctuation, (AT/T) rm3 , may be computed di- 
rectly from the two-point angular correlation function, C(9), or from the angular 
power spectrum, Ci. The former may be used if the theory provides an analyt- 
ical expression for the angular correlation function (e.g., a fitted parameterized 
expression) or, e.g., a parameterized expression for th e tra nsfer function of the 



temperature fluctuations (for transfer function see e.g 158]) 
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The expression of (AT/T) rms depends on the experimental configuration, 
e.g., beam size, switching pattern, etc. Most CMB experiments use beam- 
switching techniques, using e.g. single, double differences etc., i.e., double or 
triple beams respectively. A two-beam experiment measures the temperature 
difference between beams separated by an angle f3 on the sky 

AT/T=(T 1 -T 2 )/T, (11) 

and its expected variance is written in terms of the temperature autocorrelation 
function : 

(^L=((^)> 2,C(0, - C( » ) - <12) 

A three-beam experiment measures the difference between a field point, T/, and 
the mean value of the temperature in two directions which are separated from 
the field point by an angle j3 

AT/T = T f - i(Ti+T 2 ), (13) 

and its expected variance is given by: 



/ AT\ 2 _ // T / -1(T 1 +T 2 j 

\ T J rms \ \ 



= |C(0) - 2C(J3) + iC(20) (14) 



A four-beam experiment measures e.g.: 

AT 3 -/T = -\T S + \T j+ i - |T J+a + \T j+3 . (15) 

Where Z (Tj,Tj+i) = /3, the beamthrow, and its expected variance is expressed 
by: 

(^r) = ~ f c ^ + ~ l c ( 3 ^- ( 16 ) 

\ -l s rms 

All these expressions are true for an ideal experiment with an infinitely narrow 
beam; for a real observation it is necessary to consider the beam smearing due 
to the finite resolution of the antenna. This is usually taken into account by 
convolving the radiation intensity with a beam approximated by a Gaussian 
G{6): 

^) = ^exp(-|l), (17) 

where a is the dispersion of the Gaussian. The correlation function of the con- 
volved radiation field is given by 

C(9,a)^— £7(9V-C-/o ^ ^ (18) 



where Iq is the modified Bessel function. So for a single difference one has: 

= 2{C(0,a)-C(J3,ir)) (19) 



\ T J 

and for a double difference 



= %C{Q,a)-2C(/3,a) + i i C{2l3,a), (20) 
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where 13 is the beamthrow. For a four-beam experiment: 



/ AT \ 2 

bH =!C(0 jf r)-fC(/3 )( T) + |C(2A<T)-|C(3AcT) (21) 

\ -L / rms 

One is able to compute {^)1 ma using Equations ( |li| ) and (p^), ( p(i| ) or ( ^H ) 
according to the experimental configuration. 

The use of the angular power spectrum has in the last few years become 
the most successful representation and methodology for treatment of CMB 
anisotropics. The use of this approach is possible because the proponents of 
several models of structure formation are now making predictions of the power 
spectrum of the fluctuations. This approach uses the expression of C(6) as in 
Equation (|l^), which considers ensemble-averaged values. As an example of this 
approach consider the part icular case of a single Gaussian beam (the general case 
is addressed in Section , for comparison with the description given previously. 
In this case the smearing effect due to finite resolution of the antenna is taken into 
account in terms of a function Wi which affects the spherical harmonic expansion 
of the temperature fluctuations: 

— (a, 0,(76 ) =^2a lm Y lm {a,(j))Wi(a b ). (22) 

l ,m 

This is translated in terms of a window function Wi acting on the C; expansion 
of the correlation function: 



C(6,a) = -3- y^{2l + l)CiPi(cos9)Wi(a), (23) 

47T ' 

where [167 



47T ■ 

i 



W ; (a)= C - (2i+1) ° /8 (24) 

and 

W l {a)=e^ 2l+1 ^' A (25) 
which is usually approximated by 

Wi =exp(-/(/ + l)a 2 ), (26) 

where a = FWHM/2V2 In 2, with FWHM standing for full width at half 
maximum of the beam. Equation ^ is a particular case, for a single Gaussian 
beam, of a more general expression which is given in Section p.3| . Once given the 
Ci's one can compute the correlation using Equation (|l^) or ( |23D and AT/T from 
Equations (jT|), @ or (|l|), for the 2, 3 and 4-beam experiments respectively. 
Some experiments using switching techniques do not use a square wave chop; 
however, it may be used as an approximation. If one wants to be more precise one 
should use the appropriate beam on the sky. It is standard procedure to describe 
the details of the observation in terms of the window function, Wi, which is then 
used in the Ci expansion of (AT /T) rms in order to compute directly the expected 



temperature fluctuations: this is addressed in the following Section (Section 2. i). 
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2.3. Window functions 



A real CMB experiment is only sensitive to a given range in angular scales with 
a response to each multipole that can be characterized by the window function, 
Wi. For example, the beam smearing due to finite resolution results in the convo- 
lution of the radiation intensity with a beam approximated by a Gaussian G{9). 
In the multipole space this will correspond to a multiplication of the coefficients 
of the Legendre polynomials by the corresponding window function, as given by 
Equation (^3|) and Jig) . As already mentioned the intrinsic auto-correlations 
function d n t is transformed by convolution with the autocorrelated beam pat- 
tern of the telescope and is given by C CO nv(9) = Ci n t{0) * Ct, eam (#), where C CO nv 
is the autocorrelation function of the convolved radiation intensity, incorporating 
all the effects of the observational strategy e.g. the finite resolution of the an- 
tena, the switching pattern etc., Cbeam is the instrumental ACF defined by the 
configuration of the telescope and * refers to convolution. The beam pattern of 
the telescope depends on the experimental strategy and for example is given by 
a triple beam pattern for double differences experiments. For a real observation 
Equation (|l^) becomes 

C conv (9) = -±-J2( 2l + 1 ) ClPl ( cose ) Wl - ( 27 ) 
i 

From now on C{9) will represent the autocorrelation function of the intrinsic 
radiation intensity. The expression for the rms temperature fluctuations becomes 



/ AT\ ' 



\ J- / rms ^ n 

I 



where we are considering ensemble average values. For a Gaussian beam pattern 
the wind ow f unction describing the finite beam resolution is given by Equa- 



tion ( |2q ) 1 167] The multipole, l s , corresponding to the dispersion of the beam, a, 
is given by Z s + l/2 = l/(2sin(cr/2)) ||. This function determines the maximum 
resolution or the minimum angular scale to which the instrument is sensitive. It 
corresponds to a high I cut off which is controlled by the beamwidth of the ex- 
periment. 

The window function is given by, for a 2-beam experiment: 

W t = exp(-Z0 + IV 2 ) [2(1 - Pi(cosOS)))] . (29) 

for a 3-beam experiment: 

Wi =ex P H(Z + l)a 2 ) [±(3-4P(cos(/3)) + P(cos(2/3)))] . (30) 

and for a 4-beam experiment: 

Wi = exp(-l(l + l)o 2 ) [f - f P(cos(/3)) + f P(cos(2/3)) - |P ; (cos(3/3))] (31) 

where Bis the beamthrow. These expressions may be easily derived from Equa- 
tions (|2"i|), (^j|),(|2o|), and (pl|). The differencing techniques introduce a low I 
cutoff which is determined by the beamthrow. As mentioned above there are 
experiments that do not use a 'square wave chop', in which case the window 
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function needs to be derived separately, but in most of these cases the observing 
strategy may still be approximated as a single or double differ ence switching. For 
example, the experiments at the South- Pole (SP91 jr§, |l6§|, SP94 Q), MAX 
pl[ p4| , 186 1 and Saskatoon [211] use a sine wave chop with a weighting of ± 1 
for the South- Pole, and the first a nd s e cond har monic of the chop frequency for 
MAX and Saskatoon respectively [207, 206, 205]. All these cases may be incor- 



porated in a g eneric expression of the window function, according to White and 
Srednicki |205 [, 

The Tenerife experiment |HJ|, MSAM3 @ and OVRO p5| , |l5l[ are examples 
of a 3-beam experiment, MSAM2 J3f| is an example of a 2-beam experiment 
while PYTHON |62], 162] is an example of a 4-beam experiment. Other small 



angle experiments use different lock-in functions (ie weighting functions which 
express the weight given to each point of the beam trajectory contributing for 
a given point of the scan, for the calcu latio n of temperature , an d the overall 
normalization), e.g. South-Pole (SP91) |l6^ , [76) and ARGO [ 13^| experiments 
use a 'square-wave lock-in', the MAX e xper iment fil] , 186 1 uses a 'sine wave 
lock-in', and the Saskatoo n experiment [211] uses a double angle "cosine lock-in" 
(for details see [205, 158 1 ) . For the South-Pole experiment (SP91) the window 
function is given by: 



(32) 



where Ho(x) is the Struve function and B? is given in Equation ( pB[ ) or (|2(j). 
For the MAX experiment, considering ao = 0.65° and a — 0.425 x 0.5°, one has 
k = \/1.13 x it = 3.34, the window function is given by: 



i 

U' =k 2 J^Bt{o) V \Yi m (e,<t>)\ 2 J?(\mao\). 



21 + 1 



m — — l 



(33) 



where Ji is the Bessel function. For Saskatoon with Qo = 2.45° and a = 0.425 x 
1.44° one has k = \/l-74 x it — 4.14, the window function is given by: 



, 4-7T 

21 + 1 



2 E 



\Y lm (9,4>)\ J 2 (ma ). 



(34) 



Other experiments use interferometry techniques (e.g., CAT [171 



where the 

window function is the UV-sampling convolved with the Fourier Transform of th e 
primary beam and is usually given by the observers (see also Section ^, [ 158 1 ) . 
An easy way to parameterise the window functions is by using three numbers 
which represent the multipole at which the instrument achieves its maximum 
sensitivity, and the multipoles at which the sensitivity decreases by a factor 2 
from the peak,these are given in Table ^ of Section |^ and Table ^ of Section [| 
Shown in Fig |B| of Section are the window functions for the various observa- 
tional configurations of several experiments. According to White and Srednicki 
205, p^| , the approximation of the South- Pole experiment by a 2-beam, gives a 
difference with respect to the exact window function of the order 20%; while for 
MAX it is of the order of 10%. 



2.4- Statistical uncertainties 



The theories predict only the statistics of the temperature fluctuations and en- 
semble average values, but not their value in our universe. The C(9,a) is a 
random variable with ensemble average value given by Equation (|2^). It may be 
expressed as: 

C rm {9,a) = ^-V Qf(x)Pi(cose)Wi(cr), (35) 

l>2 



where all the quantities are as defined in Section 2.1 and [2^. The cosmic vari- 
ance, o"cos) is defined in terms of var(Qf) and represents an intrinsic limitation of 
the stati stical description of the CMB fluctuations. As already mentioned in Sec- 
tion 2.1 this intrinsic cosmic variance is different for the different angular scales 



depending on the number of independent samples in the sky on those angular 
scales. The variance is therefore larger for large angular scales. Several authors 
have considered how this uncertainty affects the estimation of the ACF. Cayon 
et al. (1991) |}3j obtained a semi-analytic expression for the probability density 
function, pdf, of the temperature autocorrelation, ACF. Scaramela & Vittorio 



(1990) [166] used the multipole expansion of the ACF and computed their dis- 
tribution by Monte Carlo simulations. Alternatively Bond & Efstathiou (1987) 
p?| , derived an analytic expression to compute a^ oa , which is expressed in 
terms of the Cv. 

var(C(0)) = JL ^(2« + l)C?/f (36) 

l>2 

Or in terms of the window function, Wi, of a given experiment: 

var(C co „„(0)) = J^(2Z + 1)C?W?. (37) 

l>2 

Equation (^) could be extended to the cases where the window function has an 
m-dependence (White et al.). Small angular scale experiments have negligible 
cosmic variance, since in this case, ergodicity applies and ensemble average is 
well represented by the angular average over the sky. The problem is that if 
the observations sample only part of the sky, the number of independent sam- 
ples in the sky at these angular scales will be reduced, this translates into an 
enha ncement of the cosmic variance, called sample variance. Scott et al (1994) 



169] derived an approximate relation between the sample variance (<r sam ) and 



the cosmic variance: 



where A is the solid angle covered by the experiment. The above methods are 
applied to experiments considering, e.g., their scan strategy, appropriate solid 
angle, etc. The sample variance may be obtained by computing the cosmic vari- 
ance via e.g. (|37]), or Monte Carlo simulations and then using the approximate 
relation given in Equation (^). We used a different approach, producing Monte 
Carlo simulations of the temperature fluctuations according to the specific char- 
acteristics of each experiment, such as the beamwidth, the sky coverage, the 
switching pattern, the beamthrow, etc. Proceeding in this way one computes 
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directly the overall theoretical uncertainty for a given experiment. This method 
was of particular interest to us because we were concerned with the confronta- 
tion of frequentist with Bayesian techniques of data analysis. These simulations 
follow different strategies according to angular scale and experimental configura- 



tion. This is described in [ 158 1 . This te chni que can also be applied to test the 



reliability of a combination of data sets [150 1 



2.5. Angular power spectrum 

Having discussed the effects of particular experimental arrangements we now 
discuss the form of the power spectrum itself. The description of the CMB 
anisotropy in terms of the angular power spectrum, Cj, has proved to be an 
invaluable method and has become a standard procedure for treatment of the 
temperature fluctuations of the CMB radiation. The competing models for the 
origin and evolution of structure, for example [ p7| , pif predict the shape and 
amplitude of the CMB power spectrum and its Fourier equivalent, the autocor- 
relation function C(ff). 

2.5.1. Power spectrum normalization 

• The galaxy clustering normalization: as 

The quantity, as, is the dispersion of the density field in a sphere of radius 
8/i _1 Mpc. The observed quantity is the variance of counts of galaxies in spheres of 
the size 8/i _1 Mpc, as(galaxies) which is related with the corresponding variance 
but in terms of mass, as (mass) via the bias parameter, b: 

b2= vl(galaxies) 

It is an observational result that as(galaxies) ~ 1 this being so, one has: b ~ 
l/as(mass). Previous to the COBE detection the most common normalization of 
models used was based on the value of as, derived from galaxy clustering assuming 



a bias b — 1 i.e. that light traces mass [104, 48j. If in fact galaxies are more highly 
clustered than matter Q , the amplitude of the initial matter perturbations 
(and hence AT rm3 /T ) necessary to produce the observed clustering is reduced 
by the factor 6. The power on large scales depends on this normalization and 
therefore depends on the Hubble constant, h, (h — Ho/100 kms -1 Mpc -1 ). It 
also depends on Qoh 2 , the quantity determining the deflection point of the power 
spectrum, where fio is the ratio of the total density of the universe to the critical 
density, for simplicty hereafter referred as to the total density of the universe. 
The relevant quantity related with this normalization is the amplitude of baryonic 
fluctuations necessary at recombination in order to obtain the fluctuations we 
see today. This amplitude lowers if the growth factor increases consequently 
lowering the normalization of the temperature fluctuations. This effect may be 
obtained, e.g., varying the value of the total density of the universe, Qo or, e.g., 
considering non-baryonic dark matter. For example the lower the value of fio 
the lower the fluctuation growth therefore the larger the normalization of AT/T. 
The non-baryonic dark matter increases the growth of fluctuations, for instance in 
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cold dark matter (CDM) models the CDM perturbations have an extra growth 
between the epoch of matter-radiation equality and the end of recombination 
oc Qoh 2 as compared to the baryonic perturbations. When baryons decouple at 
the end of recombination they fall into the potential wells of the non-baryonic 
dark matter and increase their amplitude to match the CDM perturbations. A 
flat model with non-zero cosmological constant, A, has a reduced growth factor 
when compared with a flat model with zero A. On the other hand for the same fio 
and same fit increasing the A component increments the growth of perturbations. 
There is also a dependence on the geometry of the universe via the angle-distance 
relation which affects the value of AT/T observed. This normalization is not ideal 
since it depends on the processing of the power spectrum of fluctuations and on 
the relationship between the structure we observe and the underlying density 
field. 



• The COBE-normalization 



Since the detections of CMB anisotropy other normalizations have been adopted. 
On scales > few degrees, CMB observations probe scales of 1000's of Mpc, inac- 
cessible to other means of measuring large-scale structure such as redshift surveys 
or galaxy counts. At these large angles, the structures form part of an intrinsic 
spectrum of fluctuations generated through topological defects or inflation. In 
this linear growth regime, observations of the scalar CMB fluctuations provide a 
clean measure of the normalisation of the intrinsic fluctuation power spectrum. 
This norm alisat ion has been established by a number of in dependent CMB ob- 
servations [ 17S , [To], [n]. Since the COBE detection 178] it became common 
procedure to consider the models COBE-normalised. This normalization is usu- 
ally given in terms of the quadrupole term, C2 or equivalently quoted in terms 



of the Q r 



where: 



T 



5C 2 
4tt : 



where T = 2.736°K. This quantity is obtained by fit ting 



187 



COBE data, according to several techniques j7^, [79] 
The relationship between this and the previous normalization is given by: 



(40) 

2;iven mo del t o the 
0, H etc. 



7T /Hofly (5M 



/ R 



(41) 



(Coles, 1995 where SM/M — au is the mass fluctuation observed at the 

present epoch on a scale R, and is as when R — 8/i _1 Mpc . 



2.5.2. The shape of the power spectrum 
• Doppler peaks 

Contemporary cosmological models with adiabatic fluctuations predict a sequence 
of peaks on the power spectrum which are generated by acoustic oscillations of the 
photon-baryon fluid at recombination. Photon pressure resists compression of the 
fluid by gravitational infall and sets up acoustic oscillations. The fluctuations as a 
function of k go as cos(kc 3 rj t ) at last scattering, where k is the wavenumber of the 
fluctuation, c s is the sound speed and rj, is the conformal time at recombination. 
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Figure 1. The angular power spectrum for a standard CDM model with 
ft = 1, ft b = 0.05, H = SOkms- 1 Mpc" 1 and n = 1. The high I region of 
the power spectrum exhibits the sequence of peaks generated by acoustic 
oscillations of the photon-baryon fluid at recombination while the plateau 
in the low I region is mainly due to the Sachs- Wolfe effect. At higher I we 
observe a damping of the fluctuations due to imperfect coupling of photons 
and baryons. 



This will produce a harmonic series of temperature fluctuation peaks, as shown 
in Fig. |l|, with the mth peak corresponding to k m — mn[c 3 ru. The critical scale 
is essentially the sound horizon c s rj* at last scattering [184 , |99). Of particular 
interest is the height and position of the main acoustic peak — the so called 
Doppler peak. The height depends on quantities like the baryonic content of the 
universe, Qb, and Hubble constant, Ho, whilst the position depends on the total 
density of the universe, fio and is expected to occur on an angular scale ~ 1° . The 
precise form of the Doppler peak depends on the nature of the dark matter, and 
the values of Qo, Qb and Ho. The scale l p of the main peak reflects the size of the 
hori zon at last scattering of the CMB photons and thus depends almost entirely 
]98| , 105] on the total density of the universe according to l v oc 1/VQ. This 
relation comes from the conversion of a spatial fluctuation on a distant surface 
to an anisotropy on the sky. This conversion involves the spectrum of spatial 
fluctuations, the distance to the surface of their generation and the curvature 
(or lensing in light propagation to the observer). The most important effect is 
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due to the background curvature of the universe, e.g. for an open universe a 
given scale subtends a smaller angle on the sky than in the flat universe, due 
to the fact that photons curve in their geodesies. In a ACDM universe (CDM 
model with 0) this main peak is located at around I ~ 220, approximately 
at the same scale as in a flat CDM model, with a small dependence on £1q and 



h [ 1 2 1 1 . In conventional inflationary theory [p6|, one expects the universe to be 
flat with Q = 1.0, which can be achieved if the total mass density is equivalent 
to the critical density or if there is a contribution from a cosmological constant 
A in the right amount. The height of the peak provides additional cosmological 
information since it is directly proportional to the fractional mass in baryons Qb 
and also varies according to the expansion rate of the universe as specified by the 
Hubble constant Ho; in general [Q for baryon fractions ^6^0.05, increasing Ho 
reduces the peak height whilst the converse is true at higher baryon densities. 
In the case of topological defect models, there is also a sequence of peaks on the 
power spectrum for the texture model, while for the cosmic strings model onl y 



the main Doppler peak is present with no existence of the secondary peaks [ 13C ] , 
although the situation is not entirely clear yet. 

• Damped tail 

At higher / a damping of the fluctuations occurs due to imperfect coupling of 
photons and baryons. The photons possess a mean free path in the baryons, 
A c , due to Compton scattering. Assuming a standard recombination history, the 
damping length approximately scales as \d(ti*) oc 7/2 (Qt,h 2 )~ 2 . 

• Sachs- Wolfe effect 

The low I region of the power spectrum is mainly contributed by the pot entia l 



fluctuations in the last scattering surface, the so called Sachs- Wolfe effect [ L62 ]. 
This efect corresponds to red shifting of the photon as it climbs out of the po- 
tential well on the surface of the last scattering, but also to a time dilation effect 
which makes as to see them at a different time to unperturbed photons p2| |64[ . 
Its expression is given by: 



/ATX 



1 



, , , 3^ (42) 

where S<f> is the perturbations to the graviational potential and c is the speed of 
light. For flat (il = 1) scale invariant (n = 1) models with no contribution from 
a cosmological constant this effect gives rise to a flat power spectrum, i.e., 1(1 + 
1)C; =constant, where the quantity 1(1 + 1)C; is the power of the fluctuations per 
logarithmic interval in I. Its shape is altered if one assumes a tilted initial power 
spectrum with a power law, P(k) = Ak n or if one incorporates a cosmological 
constant, or assumes other than flat models. In the last two cases another source 
of anisotropy appears, the so-called integrated Sachs- Wolfe effect and is due to the 
fact that potential fluctuations are no longer time independent. An ilustration 
of the CMB anisotropy spectrum is given in Fig. |l|. 

2.5.3. Calculation of the Ci 's Theoretical calculation of the CMB anisotropics 
are based on linear theory of cosmological perturbations. Most of these calcu- 
lations solve the Boltzmann equation using Legendre expansion of the photon 
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distribution function. This method expands each Fourier mode of temperature 
anisotropy in Legendre series up to some l ma x- This system of equations is then 
numerically evolved in time from the radiation dominated epoch until today. 



These theoretical calcul atio ns were first developed by Lifshitz (1946) [ 122 1 , while 
Peebles and Yu (1970) ]l4C| ] applied these calculations to the CMB anisotropics, 
considering only photons and baryons. Later other authors include d e.g . dark 
matter (Bond & Efstathiou 1984 1987 @, Vittorio & Silk 1984 |l9q), cur 



vature (Wilsonfe Silk 1981 ]2iq ], Sugiyama & Gouda 1992 jlSJ], White& Scott 
1995 J208[ ), tensor modes or gravitational waves (Crittenden et al. 1993 J4ql) 
and massive neutrinos (Bond & Szalay 1983 ^8j, Ma & Bertschinger 1995 |128j ], 
Dodelson, Gates & Stebbins 1995 j5S|). Other calculations have been done e.g. 
by Holtzman 1989 Q, Sugiyama 1995 pi, W ayne & Sugiyama 1995 Q which 
used an analytic approach, Stompor 1994 [p-82| 1, Seljac & Zaldarriaga 1996 172 1 
which method is based on integration over sources along the photon past light 
cone). The multipole coefficients Ci's may be obtained, for example, usin g th e so- 
lutions to the equations describing temperature fluctuations evolution by |27|. : 
Ci = k 2 dk\ At; (k, i~o)| 2 , where ATi(k,To) are the coefficients of the Leg- 

endre polynomials expansion of, AT(fc,ro), the radiation intensity fluctuations, 
t is the conformal time J dt/a(t) where a(t) is the cosmological factor, and k is 
the wavenumber of the Fourier mode. In particular in the case of large angular 
scales for fl — 1 and an initial power-law form of the density fluctuations power 
spectrum P(k) oc k n , an expression for the Ci(n) may be used: 

r = c r(f + (n-l)/2)r((9-n)/2) 
1 2 r(/ + (5-n)/2)r((3 + n)/2) 1 ' 



with n < 3, / > 2, / < 40 p7j , p4] , 165]. This expression is only adequate for 



experiments probing angular scales larger than the horizon size at recombination 
(> 1°) where only Sachs-Wolfe or isocurvature effects are dominant. In order 
to get the Ci's for smaller angular scales one can use of an alternative method; 
this consists of using the expression of the expected temperature autocorrela- 
tion function and the orthogonality pro perty of the Legendre polynomials, which 



generates Ci for all angular scales [ 15S 



In the last couple of years various groups have been producing their own CMB 
codes to generate the Ci's according to different approaches (G. Efstathiou, N. 
Sugiyama, U. Seljac and M. Zaldarriaga, M. White etc..) which are now publicly 
available. 



2.6. Likelihood analysis 

A most commonly used statistical te chni que to analyse the data is the Bayesian 
'likelihood function' analysis |l^, ^l|, which we will describe in this section. 
This analysis permits one to use the correlation information existing in the data 
and to deconvolve the switching correlations allowing one to test directly a given 
sky temperature fluctuation model. In order to apply this analysis one assumes 
a given structure formation model and considers its predicted power spectrum 
which defines the intrinsic angular correlation function: 
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Gw(0) = (^(rTi)^(r5)), (44) 

where n\, ni are two directions on the sky such that: n*i.n2=cos(#). 

This relates with the autocorrelation function of the observations by: 

/ A rpobs Arpobs \ 

G ohs {6) = (^-(»i)==^-(rii) j, (45) 

which is equal to Cint{8)*Cbeam{9) where * refers to convolution and GbeanJjf) 
is the autocorrelation function of the beam of the telescope, (see also Section 2.Z ). 
The ACF measured between two points on the sky i and j with coordinates (a, 5) 
and (a', 5') (in RA and dec), My, is defined by: 

«•-(&■)?(¥)?)■ 

These are components of a matrix (My) called the covariance matrix which 
in particular contains information about the beam switching. Adding to (My) 
the covariance matrix for the errors (-Ey) one obtains the total covariance matrix 
Vij=Mij+Eij. For Gaussian-distributed fluctuations the likelihood function for 
the total covariance matrix Vy, is 

L(A T\C int (9)) = ^= exp(-^A fr'A T). (47) 

which expresses the probability of obtaining the data set AT = 
(ATi, . . . , AT n ) given some intrinsic autocorrelation function Ci n t- 

The analysis proceeds calculating different values of dnt(O) by varying the 
parameters of the model and evaluating the respective likelihood. Searching for 
the maximum value of the likelihood function the most probable parameters of the 
model are obtained. Considering e.g. the variation of the amplitude parameter 
of the sky fluctuation model, two things can happen: the point of maximum can 
be different or equal to zero. In the first case one has a detection of fluctuations, 
in the second a monotonically decreasing function is obtained and there is not 
a detection of fluctuations in the data. The next step consists in evaluating the 
likelihood ratio L max /Lo where Lo is the likelihood evaluated at zero. This ratio 
gives an idea of how representative the maximum value of the likelihood function 
is. 

The upper and lower limits on the amplitude \JC(Q) of the fluctuations are 
obtained invoking Bayes theorem: 



P(V^(0)|AT) « L(A T\ ^/C(0))P(^/C(0)), (48) 

where P( y/C(0)\ AT) is the posterior probability (of obtaining the intrinsic sky 
fluctuations given the observed data set); L(AT\ -y/C(0)) is the likelihood (of 
obtaining the data set given the intrinsic sky fluctuations) and P(^/C(0)) is 
the prior probability. Once the form to give to the prior probability is known 
we are able to know the probability of obtaining the amplitude of the intrinsic 
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sky fluctuations given the observed data set. Some Bayesians use the Jeffreys 
prior, uniform in log(A) space i.e. Pr(A) = 1/A. This is based on the idea that 
when there is no prior information available, invariance arguments imply that 
the prior probability density of a scale factor parameter, A, (which Co is) should 
be uniform in 1/A. This poses a problem because when A — > the Pr(A) — * oo 
and the posterior probability cannot be normalized. The solution could consist 
in considering a cuto ff at low A. Ins tead we have just a ssumed a unifo rm prior 
such that: P( V /C(0)) = 1 for y/C(0) > and P( v /C(0))=0 for y'C(O ) < 0. In 
this way, the area under the likelihood curve between and \J C(0) gives the 
probability of obtaining y C(0). So looking for the value of \/C(0) such that 
this area is 95% of the total area one obtains a 95% upper limit. In a similar 
way one calculates the 95% lower limit to the intrinsic sky fluctuations (this is 
the equal tail prescription, ET) . In a similar way, the 68% upper and lower limits 
to the detected fluctuation amplitude are obtained by integrating down from the 
peak of the likelihood curve at plateaus of equal likelihood, until the area under 
the curve delimited by the upper and lower values is equal to 68 % of the total 
area under the likelihood curve (this is the highest posterior density prescription, 
HPD.) The HPD prescription, for calculation of confidence intervals, results in 
somewhat smaller upper error bars and somewhat larger lower error bars then 
does the ET prescription. Work has been done in order to develop methods to 
account for uncertainties, such as those in the beamwidth and calibration, in 
likelihood analysis of CMB anisotropy data f^j |. The beamwidth-uncertainty 
correction has a bigger effect for those data sets that probe a region of Z-space 
where the model C spectra changes rapidly with I. This is due to the fact that 
the contributions from the ±1<t beamwidth analyses no longer approximately 
compensate for each other. The calibration-uncertainty correction broadens and 
skews the probability density distribution functions toward higher values of the 
amplitude (Q™ s - ps ), and the calibration-uncertainty corrected likelihood has an 
amplitude dependent width [J74|. 



2. 7. The Gaussian auto- correlation function 

At an earlier point in CMB work, the Gaussian auto-correlation function (GACF) 
was the tool most frequently used to describe an hypothe tical sky model, and 
originally was intended to be used for small scales [ 198| . 207 1. This representation 



is not physically realistic but is useful for comparison of the results from different 
experiments. In the GACF model we assume that the intrinsic autocorrelation 
function is a Gaussian of amplitude Co and dispersion 8 C : C(8) — Co exp ( 
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even for more general autocorrelation functions is useful to define, 6 C , by 9, 

/ x 1/2 

( ~ c" (0) ) • "^Co then gives an indication of the rms amplitude of the intrinsic 
fluctuation on some coherence scale defined by 8 C . In the angular powe r spectrum 



language, this corresponds to considering C; to be expressed by [207 



Ci = 2nC 8 2 c exp(~l(l + l)6l) (49) 

The form of the autocorrelation function after convolution with the beam from 
a single instrument horn, assumed to be a Gaussian, is given by the convolution 
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of the intrinsic autocorrelation function, dt, with the beam autocorrelation 
function, C eam : C(8,a) = C„ t (#) *C(, eam (#) which becomes: 

c ^-2§Tol^ (2(2^1)- ( 5 °) 

Experiments observing the same intrinsic sky fluctuation do not, in general, 
measure the same rms signal level because this signal depends on the obser- 
vational strategy. This rms signal can only be interpreted if accompanied by 
information about the window function of the observation. In a GACF analy- 
sis we assume a power spectrum with a Gaussian shape and take into account 
the peak and area of the window function. This model is a good approximation 
for an experiment where the power spectrum is smooth and the window func- 
tion can be considered to be a Gaussian, because the two convolved spectra are 
very similar, in which case the quoted Co gives a measure of the intrinsic sky 
fluctuation. The sensitivity of an experiment to fluctuations on various scales 
is commonly described by a plot of Co vs C , assuming that the sky correlation 
function is a GACF. The experiment is most sensitive to a GACF which peaks at 
the maximum of its window function, and the amplitude of fluctuations to which 
the experiment is sensitive is parameterized by the minimum value of Co- 

In many instances experimenters now report results in terms of flat bandpower 
estimates, where the form of the angular power spectrum C; is represented by a 
flat spectrum over the width of a given experimental window. This is described 
in detail in Section la. 



3. The experiments 

This Section gives an account of the current experimental status and of the main 
sources of contamination of data from CMB anisotropy experiments. Usually 
the experiments are classified according to the angular scales they probe. The 
usual classification considers large angular scales as being angles larger than 
the horizon size at recombination i.e. 9>2°; medium scales comprehends the 
range W'<8~2° and covers the region of the Doppler peaks; small scales lie in 
the range l'<8<10' and very small scales for When analysing the data 

from CMB anisotropy experiments one must take into account contaminants like 
free-free and synchroton radiation from the Galaxy which are dominant a t low 



frequencies (~ 10-50 GHz) and dust at higher frequencies (> 100 GHz) [116]. 
These are plotted in Fig El Anot her source of contamination is the possibility 
of discrete source contamination |P_lq]. The main parameters of several recent 
experiments are displayed in Table 111. Some of the text of this secti on is base d 



upon previously published reviews by Lasenby and Hancock (see e.g [114, 102 1 ) 



3.1. Foreground emission 

3.1.1. Diffuse Galactic emission A detailed study of the diffuse Galactic emis- 
sion has been done at low frequencies (38^^1420 MHz), using radio telescopes 



11J] and at the higher frequencies of 60/^m and 100/oti using the IRAS satellite. 
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Figure 2. Estimates of the amplitudes of the anisotropic component of 
foreground emission (Anthony Lasenby private communication). As a ref- 
erence, a typical 10 [iK CMB signal is shown. 



The FIRAS and DMR instruments on board of the COBE satellite provided fur- 
ther information on the Galactic emission on 7° scales. In the microwave regime 
the emission from our Galaxy is composed of synchroton emission from cosmic 
ray electrons, free-free (bremsstrahlung) radiation from thermal electrons and 
thermal emission from dust. Fig. |^ shows estimates (Lasenby priv. comm.) of 
the relative antenna temperature contributions expected from each of these dif- 
fuse Galactic foregrounds, as a function of observing frequency v. It represents 
the anisotropic component of the emission at high Galactic latitude (|6| > 40°). 
The CMB radiation has a Planckian spectral form: 

AT A = f(v)AT = p±±- > (51) 

where x = hv/kT and AT a is the antenna temperature contribution from a black- 
body temperature AT. The spectra of the foreground emission can be given by 
the relation Ta = \ where (3 is the temperature spectral index. The different 
frequency spectral signatures of these components allow their discrimination via 
observations at several frequency ranges. The synchroton radiation is the domi- 
nant Galactic component at i^30 GHz, which results from relativistic electrons 
spiralling in the magnetic field of the Galaxy. Assuming that the number density 
of the electrons is of the form N(E) oc E j3 = (7 + 3)/2, the intensity of 
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synchroton radiation, I{v), is given by: 

I(u) oc B^V^, (52) 

where B is the magnetic flux density, and (3 is the temperature spectral index. 
The spectrum steepens from (5 ~ 2.75 at 1 GHz to /3 ~ 2.9 at 15 GHz, due 
to increased loss of energy by electrons at higher energies. A spatial variation 
of (5 is expected as a consequence of the spatial variation of the magnetic field 
||. So the modeling of this foreground involves the steepening of the spectrum 
with frequency and the spatial dependence of the frequency spectral index J3. 
The low frequency surveys can be used to extrapolate (see also Banday and 
Wolfendale 1991 Efl|) the synchroton component to other frequencies. There are 
low frequency surveys conducted at 408 MHz Q, 820 MHz Q, and 1420 MHz 



152J. The 408 MHz is an all-sky survey while the 1420 MHz covers only the 
declination range —19° < 5 < 90°. Given the frequency range probed by these 
surveys it is difficult to estimate the spectral steepening at higher frequencies. 
Bennet et al. Jh| made an attempt to account for the steepening by assuming 
that the spectrum of the local electrons reflects that of the electrons producing 
the observed synchroton emission. 

The free-free emission is also a problem at v < 50 GHz and possibly higher 
as well. Free-free emission (bremsstrahlung) results from the acceleration of elec- 
trons (when they interact with the warm ionized component of the Galaxy) in 
the electric field of the ions. Its spectral index is a weak function of frequency 
Jl6[ and is given by: 

® i} = 2 + 10.48 + 1.5 ln(T e /8000K) - In */(GHz) ' ^ 

where T e is the electron temperature and ~ 2.1 at the Tenerife experiment 
frequencies. This component of emission is large within the Galactic plane, but 
at ^<15 GHz and at high Galactic latitude it is expected to be much lower than 
the synchroton emission. It is thought that the free- free emission reaches ~ the 
same level of the synchroton at frequencies ~ 20-30 GHz, and to be larger for 
higher frequencies, although the accurate amplitude of this Galactic foreground 
is not yet well understood. 

At higher frequencies, larger than ~ 100 GHz, the contribution to foreground 
emission comes from the Galactic dust emission. The modeling of this foreground 
is difficult because it requires the knowledge of the temperatures, emissivities and 
spatial distributions of the dust grains in the Galaxy. Before the publication of 
COBE data, Banday and Wolfendale 1991 || concluded that IRAS observations 
provided useful information at high frequencies but its extrapolation for frequen- 
cies where the CMB is observable is difficult given the uncertainties in the dust 
model. Bennet et al. 1992 H, using empirical fits to the COBE FIRAS data 
estimated the dust contribution to the COBE observing frequencies and found 
that the rms dust signal was less than 8 piK at 90 GHz. Dust emission increases 
with increasing frequency with a spectral index /3 ~ —1.5. 

• Discrete source emission 
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Another source of contamination of CMB observations is the discrete source 
emission. Most of these sources have a fiat or falling spectra such that their flux 
decreases with increasing frequency. The identification of the brightest sources 
at low frequencies is made using the existing catalogues, but in some cases it 
is necessary to use separate high resolution surveys. The antenna temperature 
contribution for a source of flux, S, in a single beam is given by: 

T- - (54) 

where A e ff = n ( 2 1 /6> F wh — ) 2 ' * s ^ e en?ec tive area of a single beam. Even for a 
fiat spectrum i.e. S oc v , Ta decreases with frequency as u 2 , and is inversely 
proportional to the square of the beamwidth. On scales 6 < 0.5° at v < 20 GHz, 
the foreground radio sources are important contaminants of CMB observations. 
Franceschini et al. H concluded that for 9 < 1° and 20 GHz< v <200 GHz 



the contribution to AT/T is just below 10 5 . The Kuhr et al. catalogue [ 108 1 of 
radio sources provides information about surveys ranging in frequency from 26 
to 90 GHz. 

• Atmospheric emission 

Atmospheric emission causes serious problems to CMB observations. Water 
vapour emission lines dominate at 22 and 182 GHz, while oxygen line emission 
is significant at 60 and 118 GHz. These emission lines define the range of fre- 
quencies of the atmospheric windows through which CMB observations can be 
made. Experiences with differential observing techniques are only affected by 
the anisotropic component of the atmospheric emission. Fig. |2| shows an esti- 
mate of this anisotropic component at a good ground-based observatory. The 
atmospheric emission increases with frequency and is a serious obstacle to CMB 
observations with increasing frequencies. 



3.2. Large angular scales 

The large scale CMB anisotropy is mainly due to the so-called Sachs- Wolfe effect 
(see Section |^) . Fluctuations in scales larger then the horizon size at recombina- 
tion retain their primordial characteristics since they have not been changed by 
any causal processes inside the horizon before recombination. So the CMB power 
spectrum mirrors that of the initial seed perturbations and therefore reflects the 
primordial unprocessed power spectrum. Observations at these scales allow us to 
determine the normalization and slope of the primordial power spectrum. They 
may also be used to distinguish between the adiabatic and isocurvature fluc- 
tuations. These tasks can be complicated by the existence of a gravity wave 
background. At these scales it is not expected to obtain information about the 
Gaussianity of the fluctuations since the beam will average over the number of 
defects and the central limit theorem states that the result will be Gaussian. 
There are three main experiments operating on these scales: COBE, Tenerife, 
MIT/FIRS. 
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ACME/HEMT (G) 


UCSB 


1.5/2.1 


30/40 


MAX (B) 


Berkeley/UCSB 


0.5/1.0 


110, 180, 270, 360 


MS AM (B) 


GSFC/Chicago/. . . 


0.5/0.6 


180 + 3 higher 


White dish (G) 


CARA 


0.18/0.47 


90 


Python (G) 


CARA 


0.75/2.75 


90 


Saskatoon (G) 


Princeton 


1.5/2.45 


26-36 


ARGO (B) 


Rome/Berkeley 


0.9/1.8 


150 + 3 higher 


CAT (I) 


Cambridge 


0.25 


13.5, 15.5, 16.5 



Table 1. Some recent CMB anisotropy measurements. (S) Satellite, (G) 
Ground, (B) Balloon, (I) Interferometer, (Qbeam) Beamwidth, (Q throw) 
Beamthrow. 

3.3. Medium angular scale observations 

At these angular scales models predict larger amplitudes of the CMB fluctuations 
arising from the Doppler peaks (see Section These oscillatory features are 
very sensitive to the details and are model dependent, therefore their observation 
should allow good constraints to be established on model parameters. Comparing 
these experiments to those on larger angular scales one may separate the scalar 
fluctuations from any possible gravity wave background. At these scales the 
effect of the cosmic variance is small but a new source of uncertainty arises 
due to the small size of the sky areas probed by these experiments, the sample 
variance. According to Coulson et al. M4, at these scales, it should be possible to 
detect non-Gaussianity induced by defect models. Reionization of the Universe 
at some epoch can smear out these fluctuations to scales depending on the epoch 
of reionization. If it happened early enough it could erase perturbations up to 
scales of less than few degrees. There are several experiments operating on these 
scales SP/ACME-HEMT, Saskatoon, ARGO, Python, MAX, MS AM, etc. 

3.4- Small and very small angular scales 

At these angular scales the amplitudes of the fluctuations are expected to be 
reduced due to damping effects such as the Silk damping and damping due to the 
finite thickness of the last scattering surface (see Section 2). They may be further 
reduced by secondary scattering processes and damping mechanisms, although if 
there has been reionization the generation of CMB fluctuations at small scales is 
expected (Ostriker-Vishniac effect). These experiments are sensitive to imprints 
on the CMB from the seeds of galaxies and clusters of galaxies. Observations 
at these scales give information about the non-linear physics of galaxy formation 



21 



and about the thermal history of the universe. They may be compared with other 
experiments in order to constrain the angular power spectrum of a given model. 
Examples of experiments operating on these scales are, among others, OVRO, 
VLA, etc. 

3.5. Future Experiments 

In the near future it is expected that new instruments will be built with better 
accuracy and a broader range of frequency coverage, providing improved quality 
CMB data. For balloon experiments one can improve the quality of the data 
through the use of long-duration balloon flights, e.g. launched in Antartica, 
and circling the Pole, and pixels sampled, and the use of arrays of detectors in 
order to extend the frequency range observed. Examples of these are Boomerang, 
Maxima, and TOPHAT. On the ground a great improvement is expected from 
interferometers, one example is the Very Small Array (VSA), to be built by 
Cambridge and Jodrell Bank in the U.K., and to be sited in Tenerife, which 
is expected to be operational by the year 2000. The CAT is a prototype for 
this more advanced instrument. The VSA is expected to give detailed maps 
of the CMB anisotropy with a sensitivity ~ 5 /iK and comprehending a range 
of angular scales from 10' to 2°, and covering a frequency range of 28 and 38 
GHz. This instrument uses an interchangeable T-shaped configuration of 10- 
15 horn elements and simulations have shown that such a configuration is able 
of attaining a good sensitivity over the range of angular scales planned to be 
used. It is also expected to measure the values of Q and Ho with an accuracy of 
better than 10% due to the good accuracy expected over the region of the first 
and secondary Doppler peaks. Simulations have shown that this instrument will 
be quite sensitive to non-Gaussian features on these angular scales. For future 
satellite experiments one has the Microwave Anisotropy Probe (MAP) and Planck 
Surveyor satellite. MAP has been selected by NASA as a Midex mission and is 
expected to be launched between 1999 and 2001 while Planck has been selected 
by ESA as an M3 mission and will be launched probably around 2005. A very 
important characteristic of satellite experiments is that they are not affected by 
the problems caused by the atmosphere. Consequently a satellite is capable of 
full-sky coverage and has the potential to map features on large angular scales 
(>~ 10°). On the other hand it has more problems in reaching resolution on 
smaller angular scales due to the limitations imposed on the dish size. The MAP 
median resolution of its channels is around 30 arcmin while the best angular 
resolution is 18 arcmin, in its frequency channel at 90 GHz. Consequently it may 
have problems in determining the shape of the first and almost certainly of the 
secondary Doppler peaks of the angular power spectrum. Planck is expected to 
attain a resolution around 4 arcmin, with a median resolution of its 6 channels 
of about 10 arcmin. Consequently it will be possible to determine the angular 
power spectrum with good accuracy including the secondary peaks, and therefore 
the determination of the cosmological parameters with good accuracy. With 
the good angular resolution of Planck surveyor it is expected to obtain a joint 
determination of and Ho to 1% accuracy. All these, balloon, ground-based 
and satellites experiments to come constitute a good improvement of CMB data 
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and represents an important step towards understanding the characteristics of 
our Universe. 



4. The Tenerife experiment 

We here give a summary of the analysis and interpretation of results obtained on 
relatively large angular scales of ~ 5° by the Tenerife experiment. Observations of 
fluctuations in the Cosmic Microwave Background (CMB) have been widely rec- 
ognized to be of fundamental significance to cosmology, offering a unique insight 
into the physical conditions in the early Universe. As we have seen the ampli- 
tudes and distribution of such fluctuations provide critical tests of the origin of 
the initial perturbations from which the structures seen today have formed. On 
scales > few degrees, CMB observations probe scales of 1000's of Mpc, inaccessi- 
ble to conventional astronomy. At these large angles, the structures form part of 
an intrinsic spectrum of fluctuations generated through topological defects or in- 
flation. In this linear growth regime, observations of the scalar CMB fluctuations 
provide a clean measure of the normalisation of the intrinsic fluctuation power 
spectrum. This nor mal isation has been established by a number of independent 
CMB observations Jl78| , fro| , |9l| . In many theories, tensor CMB fluctuations from 
a background of gravitational waves are also expected to be significant on these 
large scales, and measuring the slope of the power spectr um offers the potential 
to constrain this contribution to the CMB anisotropy (91, 18C , [IB]. A comparison 



of the large-scale anisotropy results with those on medium scales can also be used 
to separate the scalar and tensor components under the assumption of a specific 
cosmological model. 

The Tenerife CMB experiments were initiated in 1984, with the installation of 
the first 10 GHz switched-beam radiometer system at the Teide Observatory on 
Tenerife Island. A subsequent programme of development has led to the present 
trio of independent instruments working at 10, 15 and 33 GHz. The ultimate 
objective is to obtain three fully sampled sky maps covering some ~ 5000 square 
degrees of the sky at each frequency and attaining a sensitivity of ~ 50 /iK at 
10 GHz, and ~ 20 fiK in the two highest frequency channels. Drift scan ob- 
servations have been conducted over a number of years covering the sky area 
between Dec=+30° and +45°. The deepest integrations have been conducted in 
the Dec=+40° region and resulted in str ong evidence for the presence of individ- 
ual CMB features (Hancock et al. 1994 

Davies et al. 1995 [|5l| (hereafter Paper I) described the performance of the 
experiments and gave an assessment of the atmospheric and foreground contri- 
butions to our data at Dec=+40°; Hancock et al. 1996 j8s| (hereafter Paper II) 
analysed in detail the results and cosmological implications of such observations. 
Here we present a sum mar y of the analysis and results published in Hancock 



et al. 1996 J8£ ]. Section 4.1 gives a brief description of the observational strategy. 
In Section VL2 we use one of the several statistical methods to calculate the level 
of the detected signals and their origin. A statistical comparison with the results 
of the COBE DMR two-year data conducted by Hancock S. and Tegmark M. is 
mentioned in Section VLA and used to place limits on the spectral index of the 
primordial fluctuations. 
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4-1. The scans at Dec +40° 



4-1-1 . Observations Observations were conducted at the three frequencies 10, 
15 and 33 GHz by drift scanning in right ascension at a fixed declination of 40°. 
The measurements were made independently at each frequency, using separate 
dual-beam radiometer systems as described in Paper I. The three instruments 
are physically scaled so as to produce approximately the same beam pattern 
(FWHM~ 5°) on the sky, thus allowing a direct comparison of structure between 
frequencies. A characteristic triple beam profile (switching angle 8°1) is obtained 
by the combination of fast switching (63 Hz) of the horns between two indepen- 
dent receivers and secondary switching (0.125 Hz) provided by a wagging mirror. 
We make repeated observations of the sky, binning the data in 1° intervals in 
RA and stacking them together in order to reduce the noise as compared with 
individual measurements. As a consequence of using two independent channels, 
the receiver noise contribution to the final data scans is reduced by a factor v2 
compared with single channel observations. The data considered here are those 
presented in the preliminary report by Hancock et al. (1994) J9l[] and Hancock 
et al. (1996) Q. We restrict our analysis to the RA range 161° - 230° cor- 
responding to Galactic latitude b > 56°. At these high latitudes foreground 
emission from the Galaxy is expected to be at a minimum. This sky region has 
also been selected (see Paper I) to be free from discrete radio sources above the 
1.5 Jy level at 10 GHz. The contribution of unresolved radio sources is expected 
to be AT/T < 5 x 10~ 6 at 15 GHz and significantly smaller at 33 GHz |, [57). 
The revised sensitivities per beam-sized area are 61, 32, 25 and 20 fiK at 10, 15, 
33 and 15+33 respectively. For details see Paper II. 



4-1-2. Reliability of the detected signals Determination of the amplitude of 
the CMB component of the structure requires one to consider the contributions 
of random noise and foreground signals to the observed data scans. The former 
has its origin in the thermal variations in the receivers and in the fluctuating 
component of the atmosphere, whilst the latter consists primarily of free-free and 
synchrotron emission in the Galaxy, plus emission from the Sun and Moon. Of 
these effects, only the Galactic emission remains constant from day to day at a 
given frequency. In Paper I it was estimated a maximum Galactic contribution 
of AT rms = 4/iK in the results at 33 GHz. An improved separation between the 
Galactic and the cosmological signal at each frequency has been under consider- 
ation via a new maximum entropy based meth od to reconstruct the intrinsic sky 
fluctuations as developed in Jones et al. 101 1 . 

The presence of correlated atmospheric component was taken into consider- 
ation in paper II, via a new split of the 33 GHz data into two subsets X and Y 
such that both channels of a given scan are included in the same data subset. An 
analysis similar to that presented in Hancock et al. (1994) [^l[ gives an astro- 
nomical signal (pi = &fx+Y)/2 ~ a \x-Y)j'i) with an amplitude a s = 43 ± 12 uK. 
The value of the signal quoted in Hancock et al. |)lj was a a id — 49 ± 10 fj,K; the 
difference with the improved estimation is certainly due to the subtraction of the 
atmospheric signal. The difference between both estimates (a = (o1 ld — al) 1 / 2 ) 
is ~ 23 uK which is our best assessment of the atmospheric contamination in the 
analysis based on the addition and difference of the 33 GHz data; this value is 
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in good agreement with estimates obtained using other methods. For details see 
Paper II. 



4-2. Statistical analysis 

The estimates of the astronomical signals made in the previous section can be 
improved on using a detailed statistical analysis. Here we use a likelihood analysis 
and take into account the contribution of the correlated atmospheric noise by 
enlarging the error bars on the stacked scans. These modified scans are also used 
to study the presence of common features at 15 and 33 GHz by the calculation 
of the cross-correlation function (for details see Paper II). Here we will focus our 
attention to the likelihood analysis technique. 



4-2.1. Likelihood analysis This analysis takes into account all the relevant pa- 
rameters of the observations: experimental configuration, sampling, binning, etc. 
(see Section |^). The combination of the atmospheric and instrumental noise can 
be modeled by a Gaussian distribution uncorrelated from point to point, (Pa- 
per II), which implies that the noise only contributes to the diagonal terms of 
the covariance matrix. We also assume that the astronomical signal is described 
by a Gaussian random field and therefore our results correspond to a superpo- 
sition of Gaussian fields in which all their statistical properties are specified by 
the covariance matrix, which takes into account the full correlation between the 
data points. The likelihood analysis procedure is to vary the parameters of the 
model to obtain different values of C; n t(#) and then to calculate the likelihood 
of observing the data set AT. Then the most probable model parameters are 
determined and a detection of fluctuations will show itself as a peak in likelihood 
away from zero (see Section ^ . 

• Fluctuations with a Gaussian auto-correlation function (GACF) 

We have analyzed our results for two hypothetical sky models, the first of 
which corresponds to a signal described by a Gaussian auto-correlation function 
(ACF) with amplitude y/Co and width 8 C (see Section ^ ). This i s no t a realistic 



physical scenario but has been used widely in the past M, 151, 198] because it 



provides for an easy comparison between the results of experiments with different 
configurations. The intrinsic ACF for these models is a Gaussian of amplitude 
\ Z(7o and dispersion 9 C , which is modified accordingly by our triple beam filtering 
([198 1, see Section^. Our instrument is sensitive over a range of coherence angles 



l°;S# c ;Sl0 o , attaining peak sensitivity for a coherence angle of 4°. 

We ha ve a pplied this analysis to the latest data at declination 40° described 



in Section 4.1 and Paper II. We have analyzed the X and Y subsets for 15 and 
33 GHz, the total stacked scans at these two frequencies, and our best scan 
15+33. The results for 9 C — 4° are presented in the third column of Table ^. 
The amplitude of the intrinsic signal corresponding to the maximum likelihood 
is given, along with the one-sigma confidence bounds calculated in a Bayesian 
sense with uniform prior. All results look consistent with clear detections at the 
two to three sigma level, and mean values of the signal slightly smaller than those 
presented in Hancock et al. |9]]], due to the improved estimate of the error bars 
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in the stacked scans (see Section ^l] and Paper I). Fig. |^ presents the contours of 
equal probability for the 15+33 scan. We see a well defined point of maximum 
likelihood at 8 C ~ 4° and ^/Co ~ 50 /xK. 

• Fluctuations with a power law spectrum 

The second model considered here is more interesting from a cosmological 
viewpoint. It corresponds to the prediction of the power law form (P(k) oc k n ) 
for the spectrum of the primordial fluctuations (see Section ^). Consider- 
ing only the Sachs- Wolfe part of the spectrum of the fluctuations the intrin- 
sic ACF can be expressed as C{9) = j^Y^A 21 + l)C" i (S) Pi(cos ff) with C[ S) = 
r['+(5-n)/2 j r[(3+")/2 ] wnere the sum is extended to the multipoles I < 60 
which corresponds to the range of angular sensitivity of our experiments. For 
I 20, standard models predict additional contributions to the CMB anisotropy, 
as one moves into the low I tail of the CMB Doppler peak. Hence fitting for 
the Sachs- Wolfe term alone to CMB data on these scales can lead to the derived 
values for n being increased by as much as 10% over the true primordial value. 
This point should be borne in mind when comparing the limits on n from the 



Sachs- Wolfe term (Section 13) to those from a fit to a full CDM type functional 
form as in Section M. For a given value of the spectral index n, the intrinsic ACF 
is a function only of Qrms-pa ■ 

This analysis was applied to the recent data described in Section ^jj. Fig. ^| 
shows the likelihood surface as a function of Q r ms- P s and n for the 15+33 scan. 
The peak likelihood forms a ridge displaced from zero in Q rm3 _ ps and corresponds 
to a 3 — 4 sigma detection of structure for each value of n considered. The shape 
of the surface implies that all values of n in this range are equally likely. This 
is predominantly a consequence of our observing technique which samples only a 
small angular range of the spectrum of fluctuations. Thus whilst our observations 
provide a good measure of the fluctuation power on ~ 4° scales, they do not in 
themselves contain sufficient information about the distribution of power with 
angular scale to allow a useful determination of the spectral slope: for this one 



must compare with experiments on other angular scales (see Section i.c and 
Sections |H| and ^j). For the specific case of a Harrison-Zel'dovich spectrum (n = 1) 
the results of the likelihood analysis are given in column two of Table H. The 
normalisation Q rm s- V s corresponds to the maximum of the likelihood function 
and the confidence intervals are at 68 %, calculated in the standard Bayesian 
manner using uniform prior. We see that in general the results are consistent 
and agree with a global normalization of the quadrupole Q r ms- P s ~ 20 — 25 
fj,K. Our best estimate for Q r ms- P s of 22lg°/iK from the 15+33 scan is reduced 
over the value of 26 ± 6/jK previously reported due to our now having properly 
accounted for the correlated atmospheric noise. 



4-3. Statistical comparison with CO BE DMR 

A comparison between the results of different CMB experiments offers the op- 
portunity to check independent measurements, to extend the range in frequency 
and angular scale, to constrain cosmological models and, if the sensitivity of the 
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Figure 3. Top Figure: Contour levels of equal likelihood for the 15+33 
scan in the case of a Gaussian shaped ACF. The contours correspond to 
10, 20, 30, 90 % of the total probability distribution. The Tenerife con- 
figuration obtains maximum sensitivity for coherence angles in the range 
2°<#<6°; structure is clearly detected at ~ 50^K over this angular scale 
range. Bottom Figure: The two-dimensional normalised likelihood surface 
as a function of the spectral index n and the normalisation Q rm s-ps for 
the Tenerife data (in the case of fluc2Tiations with a power law spectrum). 
The projected contours are at 68 %, 95 % and 99 % confidence. 



Table 2. Results of the likelihood analysis for a Harrison-Zel'dovich spec- 
trum of fluctuations (second column) and for a Gaussian ACF (third col- 
umn) . 
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experiments is sufficient, to compare features. There are seve ral e xperiments op- 
erating on angular scales of a few degrees: MIT @, COBE [jl7|, (n}), RELIKT 
183 ], ARGO 1 132] and Tenerife. The first comparison between independent CMB 
observations was made by Ganga et al. (1993) [^o| who found a clear correlation 
between the results of the first year of COBE DMR observations and those of the 
MIT experiment. De Bernardis et al. (1994b) ^] have also made a statistical 
comparison between the amplitude of the signal reported for ARGO and that of 
the COBE DMR first year results, from which they constrain the spectral index 
to be 0.5 < n < 1.2 in the absence of any gravity wave background. In the pre- 
liminary report (Hancock et al. 1994 |9l[|) the amplitude of the signal detected 
on ~ 5° scales in our Dec=+40° data was compared with that on ~ 7° scales for 
the first year of COBE data. It was found that both results were consistent with 
inflationary models (n>0.9) but with a favoured spectral index of n — 1.6 (this 
value decreases if we use in the comparison the new astronom ical signal corrected 
for the atmospheric contamination). Lineweaver et al. (1995) Jl25} | have presented 
the first direct comparison of CMB features between the two-year COBE DMR 
data and the Tenerife Dec=+40° observations, confirming the agreement in the 
level of the normalization of both experiments and providing clear evidence for 
the presence of common hot and cold spots in both data sets. Here we give a 
brief account of a new comparison carried out by S. Hancock and M. Tegmark 
|§9[ . This comparison is different to that conducted in Hancock et al. (1994) |}lj 
in that they use more rigorous comparison technique that utilises the likelihood 
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function to incorporate fully the effects of cosmic and sample variance, random 
noise and the interdependence of the model parameters. In addition to the im- 
provements from this revised analysis, the new results also reflect the increased 
sensitivity of the COBE data after two years of observing, along with the more 
accurate estimate of the cosmological signal in the Tenerife data. 



4-3.1. Properties of the two data sets The instrumental profile of the COBE 
data is described approximately by a Gaussian beam with FWHM~ 7°. As can 
be seen in Fig. 1 of Watson et al. (1992) there is a range of angular scales to 
which the COBE DMR and Tenerife experiments are both sensitive. Measure- 
ments taken by COBE cover the full sky, but to determine the CMB fluctuations 
the region of the Galactic plane has been excluded (\b\ < 20°) thereby introduc- 
ing a degree of uncertainty in estimating th e pr operties of the global field; this 
effect is commonly termed sample variance 169 1 (Section 2). The uncertainties 
in the COBE two year results are dominated by the effect of cosmic variance i.e. 
the fact that our stochastic theory describes the Universe as a particular realisa- 
tion of a random field (Section 2). Together the cosmic and sample uncertainties 
form an intrinsic limitation of the COBE experiment since unlike random errors 
they are not reduced by increased integration time. The two year COB E da ta 
have been analysed independently by a number of authors (e.g. O , |l7[ f78| , 2 1 2 1 ) . 
All find evidence for statistically significant structure at an amplitude consistent 
with that of 30 ± 5 /jK rms on a 10° scale announced by Smoot et al. (1992) for 
the first year data. The best fit values for n and Q r ms- V s depend on the pre- 
cise analysis techniques employed, but are generally consistent with the values of 
n = 1.10 ± 0.29, Qrms-ps = 20.3 ± 4.6 found by Tegmark and Bunn (1995) 
187] for the combined 53 and 90 GHz data with the quadrupole included. In the 
case of the Tenerife experiment the double-switching scheme removes the contri- 
bution of low order multipoles decreasing the cosmic variance of the signal on 
these large scales; the major source of uncertainty is produced by the partial sky 
coverage (sample variance) and the instrumental noise. The region observed by 
the Tenerife experiments covers ~ 5000 square degrees but here we have limited 
our analysis to our region of deepest integration at high Galactic latitude which 
constitutes a sample ~ 500 square degrees. For such a region the uncertainties 
due to the part ial sky coverage dominate over the intrinsic variance by a factor 
~ 10 [169, 158 1 and the combined uncertainty is approximately of the order of 
that introduced by the instrumental noise in the 15+33 scan. 

What is required is a data analysis technique that allows the joint proba- 
bility of any combination of the model parameters to be calculated and which 
implicitly takes into account random errors and cosmic and sample uncertaint ies. 



The Bayesian approach using the likelihood function as described in Section 4.2 



attempts to do precisely this. The likelihood function peaks at the most likely 
parameters (the best estimate of the true values if the likelihood function is unbi- 
ased) and has some distribution which through Bayes theorem is representative of 
the combined effects of the cosmic, sample and random uncertainties. The issue 
of how well this distribution reflects the true uncertainties is addressed by com- 
parison of the Bayesian probability distrib utio n with that obtained from direct 
Monte Carlo simulations of the data (see [158]). The Bayesian and frequentist 
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approaches are found to be consistent for the Tenerife data and to a good approx- 
imation the likelihood function is also seen to be an unbiased estimator of the 
model parameters. Consequently the likelihood surface for the joint Tenerife and 
COBE data set provides the definitive means of comparison of the observations 
under some assumed sky model. 



4-3.2. The Tenerife-COBE likelihood function S. Hancock and M. Tegmark 
applied the likelihood analysis to the COBE two year data and the Tenerife 
15+33 scan, assuming a power law model with free parameters n and Q rms _ ps . 
The COBE Galaxy-cut two-year map consists of 4038 pixels, whilst the Tenerife 
Galaxy cut (RA 161° — 230°) scan contains 70 pixels, requiring a 4108 x 4108 
covariance matrix for a joint likelihood analysis of the data. For details about 
the implementation of the joint Tenerife-COBE likelihood function see Hancock 
et al. 1996 §§. 

For our Tenerife data, the best estimate of the cosmological signal is obtained 
from the 15+33 combined scan after correction of the error bars for the correlated 
atmospheric noise term. The possible contribution of Galactic signals has been 
estimated from the 10 GHz data to be less than 4 fiK at 33 GHz and has not 
been considered in the current comparison. The normalised likelihood function 
for this scan, as plotted in Fig. ^, represents the joint probability of obtaining 
a given combination of n and Q rms _ ps . On its own, the Tenerife configuration 
provides less leverage on the slope of the spectrum than the COBE satellite. This 
is because the Tenerife experiment is insensitive to the largest angular scales, 
and because the one- dimensional shape of the dec+40° strip makes it difficult 
to separate the power contributions from different scales. In other words, a 
narrow strip corresponds to wide window functions in £-space, with considerable 
aliasing of small-scale power onto larger scales (just as the case is with one- 
dimensional "pencil beam" galaxy surveys). As a result, the Tenerife data can 
be equally well fit by a range of n and Q rm s- V s values, resulting in a likelihood 
ridge in (n, Q rms „ ps )-space with minimal discriminatory power for the parameter 
n. In contrast, the COBE observations are sensitive to the slope to the extent 
that the likelihood surface is peaked in the n-dimension. Combining the COBE 
information with the Tenerife data improves the situation in two ways: it extends 
the lever arm on the spectral slope from the COBE scales down to the 4° scale of 
Tenerife, and in addition eliminates the above-mentioned aliasing problem, since 
the joint data set is no longer one-dimensional. 

Fig. [l| shows the confidence contours obtained from Bayesian integration un- 
der the combined COBE and Tenerife likelihood surface assuming a uniform 
prior. The 68% joint confidence region in (n, Q rms „ ps )-space encloses a region 
from 0.90 to 1.73 in n for Q rms _ ps in the range from 12.1 to 22.9 /iK, with 
the peak at n = 1.37, Q r ms- V s = 16.1/iK. Marginalizing over Q rm s~ V s with a 
uniform prior, one obtains the probability distribution for n as given in Fig. 
corresponding to n = 1.33 ± 0.30 at 68% confidence. The resulting limits on the 
normalization, conditioned on n = 1 as is customary, are Q r ms - ps = 21.0 ± 1.6. 
The corresponding results in Tegmark and Bunn (1995) 187], hereafter "TB95" 
using just the COBE data, and including the weak correlated noise term 
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were n = 1.10 ± 0.29 and Q rm s- P s = 20.3 ± 1.5, with the peak likelihood located 
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Figure 4. Top Figure: Constraints on the quadrupole Qrms-ps and on the 
spectral index n of fluctuations obtained from a joint likelihood analysis 
of the Tenerife data and COBE DMR two-year data. The contour levels 
represent 68 %, 95 % and 99 % of the region of joint probability. The peak 
of the distribution lies at n — 1.37, Q r ms~ ps = 16.1^iK and is identified by 
the cross. Bottom Figure: The marginal likelihood for the spectral index 
n as obtained from the joint analysis of the Tenerife and COBE data. The 
spectral index is seen to lie in the range 1.02 < n < 1.62 at 68% confidence, 
with a best fit value of n — 1.33. 
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at n = 1.15, Qrms-ps = 18.2^tK. In other words, although the total normaliza- 
tion has risen by a mere 3%, the slope estimate has risen by 12% and the peak 
likelihood has been shifted to higher n and lower Q rm s-ps- This indicates that 
the higher angular resolution data from the Tenerife experiment contains slightly 
more power on small scales. As explained in Section E^q, this is not unexpected, 
since the presence of a Doppler peak would cause a rise in the power spectrum 
at higher I. 

The large angular scale CMB observations from Tenerife and COBE probe 
fluctuations that have yet to go non-linear and the shape of the power spectrum 
is thus insensitive to the exact abundance of the baryonic mass (Qb) an d to the 
value of h = Ho/ 100 km s -1 Mpc -1 . However, together the Tenerife and COBE 
observations provide a direct measure of the CMB power spectrum normalisa- 
tion, against which one can compare intermediate scale observations and hence 
discriminate between competing cosmological models. Here we stress the impor- 
tance in the agreement between the derived normalisation for the independent 
Tenerife and COBE experiments, which are subject to different systematic errors 
and different foreground contamination. The accuracy to which we know the 
normalisation of the power spectrum clearly becomes an issue when comparing 
with smaller scale observations to determine cosm ologi cal parameters. This is 
a particular concern if, as has been suggested [ fls} |l8Cl| , |l[, a component of the 
large scale anisotropy may be due to tensor metric perturbations produced from 
a background of gravitational waves. These can arise naturally in inflationary 
scenarios and would contribute a component C, =< la^rj 2 > to the observed 

(T) 

CMB angular power spectrum. The ratio of the tensor modes C, to the scalar 
modes C; (primordial density fluctuations) is highly suppressed for fluctuations 
contained within the horizon volume at recombination and hence their contribu- 
tion is only significant on scales > 2°. Consequently the existence of a tensor 
contribution has implications when comparing the large scale anisotropy level 
with that on smaller scales and with large scale structure observations in order 
to test cosmological models. The anisotropy measurements from Tenerife and 
COBE fix the sum of C} T 'and C\ , but the separation of the two terms requires 
a comparison with smaller scale o bserv ations under the assumption that a given 
cosmological model is correct Jf5| , 180 . In the case of power law inflation a re- 
lation exists between the tensor to scalar ratio C% I G2 anc ^ the s l°P e of the 



primordial power spectrum {l5| , 18C ] 



Cf>/Cj 5) «7(l-n) (55) 

from which we see that the two contributions are comparable at n = 0.85 with 
Cy decreasing relative to C\ for higher values of n. Thus the limit of n > 1.0 
obta ined from the analysis of the combined Tenerife and COBE data in Sec- 
tion 1.3.2] implies that it is unlikely that the tensor component will be dominant 



in such models. This is investigated in more detail by comparing with medium- 
scale anisotropy results. 

The work required to place the existing medium scale anisotropy results into 
a common statistical framework and then to comgare them with the predictions 
of cosmological models, is carried out in Sections 



pare men 
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4 -4- Conclusions 



We have presented a summary analysis of the measurements taken at 
Dec=+40° by the Tenerife experiments at 10, 15 and 33 GHz. After accounting 
for both local atmospheric and discrete radio source foregrounds, the 15 and 33 
GHz data at high Galactic latitude are seen to contain statistically significant 
signals which have their origin in common hot and cold features. The cross- 
correlation function between the data at these two frequencies demonstrates that 
the amplitudes and shapes of the structures detected at 15 and 33 GHz are sim- 
ilar (paper II). This, combined with our measurements at the lower frequency of 
10 GHz, implies that the CMB signal dominates over the Galactic contribution 
at 15 GHz, and that the maximum possible Galactic contribution at 33 GHz is 
smaller than 10 % of the detected signal. Our best estimate of the cosmological 
signal is Q rms -ps = 22^g°/iK for an n — 1 inflationary spectrum. This ampli- 
tude is reduced by 4^iK over that previously reported for the same data set and 
results from an improved separation of signal from atmospheric noise. Compar- 
ison of the Tenerife and COBE two year anisotropy detections by means of the 
likelihood function allows a detailed investigation of the allowed parameter space 
for a power law model of the fluctuation spectrum. The best fit values of n and 
Qrms-ps are 1.37 and 16/iK, and marginalising over Q r ms- P s we find both data 
sets consistent with n in the range 1.0 < n < 1.6. These results support infla- 
tionary models, which predict n ~ 1 and future improvements in the Tenerife 
data j88[ and use of the 4 year COBE data will narrow the range of allowed n. 
Improvements of this kind are important to determine the power spectrum nor- 
malisation, since only on these large angular scales is it possible to place limits on 
the tensor to scalar ratio independent of the precise details of the cosmological 
model. In Sections |E]^3] improved limits have been placed on n, by combining the 
large scale anisotropy measurements from Tenerife and COBE with observations 
of medium scale anisotropy measurements, but at the expense of assuming an 
underlying cosmological model (in our case CDM). 



5. Constraints on the Cosmological Parameters using current 
CMB observations-I 

Observations of the Cosmic Microwave Background (CMB) radiation provide in- 
formation about epochs and physical scales that are inaccessible to conventional 
astronomy. In contrast to traditional methods of determining cosmologica l pa - 



rameters, which rely on the combination of results from local observations [ 13£ ] , 
CMB observations provide direct measurements | p7 , 207 over cosmological scales 



thereby avoiding the systematic uncertainties and biases associated with conven- 
tional techniques. As mentioned in Section 2 the principal cosmological informa- 



tion is contained in the acoustic peaks [ |98| , 170 1 in the power spectrum, which 
are generated during acoustic oscillations of the photon-baryon fluid at recombi- 
nation Jm)]. The main acoustic peak, known as the 'Doppler peak', is a strong 
prediction of contemporary cosmological models with adiabatic fluctuations and 
is expected to occur on an angular scale ~ 1° (see Section 2); its observation is 
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a major goal of observational cosmology. In the case that this peak was not ob- 
served, this would imply either that medium-scale primordial CMB fluctuations 
had been wiped out by reionization |m| or that there is a fundamental flaw in our 
theory. On the contrary, a conclusive observation of the first peak would provide 
strong support for current theoretical models and the determination of its angu- 
lar position would constitute a direct probe of the large scale geometry of the 
Universe. The angular scale l p of the main peak reflects the size of the ho rizon 
at last scattering of the CMB photons and thus depends almost entirely |)8| , 105 1 



on the total density of the Universe according to l p oc l/VTl. In conventional 
inflationary theory flssfl, one expects the Universe to be flat with S7 = 1.0, which 
can be achieved if the total mass density is equivalent to the critical density or if 
there is a contribution from a cosmological constant A. The height of the peak 
provides additional cosmological information since it is directly proportional to 
the fractional mass in baryons fij, and also varies according to the expansion rate 
of the Universe as specified by the Hubble constant Ho; in general bq] for baryon 
fractions f2(,;S0. 05, increasing Ho reduces the peak height whilst the converse is 
true at higher baryon densities. Furthermore, by measuring the amplitude of the 
intermediate scale CMB fluctuations relative to those on large scales it is possible 
to place tight limits on the spectral slope n of the initial primordial spectrum of 
fluctuations. The latter is predicted by inflationary theory to be approximately 
scale invariant, in which case n ~ 1 .0, although the presence of a background 



of primordial gravity waves 18C] would generally lead to lower values of n. 



Thus, in summary, by comparing large and intermediate scale CMB observations 
and tracing out the Doppler peak, it is possible to directly estimate Q, fib and Ho 
and to probe inflationary theory and the existence of primordial gravity waves. 
Recent improvements in the quality of CMB data, in particular on the angular 
scales probed by the CAT and Saskatoon experiments, now make this exercise of 
great interest. 

Here we pre sent the analysis and results published in Hancock et al. 1997 
H. In Section 



5.1 



we describe the method applied for the intercomparison of 



experiments and models. In Section 5.2 we proceed to the discussion of the results 



obtained. For more information see also 
5.1. Method 
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Clear detections of CMB anisotropy ha ve b een reported by a number of different 
groups, including the COBE satellite | 178, E"i |; ground-based switching exper- 



iments such as Tenerife plf , Python [ 162 1 South Pole B] and Saskatoon 



1391; balloon mounted instruments such as ARGO {52) , MAX 118611 , and MS AM 
35L B6J and more recently the ground-based interferometer CAT |17l| (these were 
the most recent detections at the time of publication of paper Hancock et al. 1997 
( @1 ) ) • Recently other new detections have been reported and are considered and 
included in the analysis carried out in Section ^[ Given the difficulties inherent 
in observing CMB anisotropy, it is possible that some of these results are con- 
taminated by foreground effects and it is clear that determining the form of the 
CMB power spectrum in order to trace out the Doppler peak requires a careful, 
in-depth consideration of the CMB measurements from the different experiments 
within a common framework. 
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We consider all of the latest CMB measurements at the time of the analysis 
described in paper Hancock et al. (|h]]), including results from COBE, Tenerife, 
MAX, Saskatoon and CAT, with the exception of the MSAM results (see below) 
and the MAX detection in the Mu Pegasi region which is contaminated by dust 
emission |3f|. (In Section ^ this data set is updated by inclusion of some more 
recent detections and a comparison of the results is there given.) The competing 
models for the origin and evolution of structure predict |27], K)8|, the shape and 
amplitude of the CMB power spectrum and its Fourier equivalent, the autocor- 
relation function C(8) =< AT(ni)AT(ri2) > where ni ■ 112 = cos9. According 
to Section |^ expanding the intrinsic angular correlation function C(8) in terms 
of spherical harmonics one obtains C(8) = X/i>2^' 1) C 1 Pi (cos 8) /4ir, where 
low order multipoles I correspond to large angular scales 8 and large Z-modes 
are equivalent to small angles on the sky. The Ci's are predicted by the cosmo- 
logical theories and contain all of the relevant statistical information for models 
described by Gaussian random fields |27| • As mentioned in Section ^| the different 
exp eriments s ample different angular scales according to their window functions 



Wi [206, 205]. The window function Wi specifies the relative sensitivity of an 
experiment to a given i-mode, and the observed power in CMB fluctuations as 
seen through a window Wi is given by 

A 2 00 

C ohs {Q) = = 5^(21 + l)CiWi/4rr. (56) 

i>2 

Given Wi, then for the Ci's corresponding to the theoretical model under con- 
sideration it is possible to obtain the value of AT„t s one would expect to observe 
using the chosen experiment (see also Section ^). This value can then be com- 
pared to the value actually observed to test the cosmological model. Shown in 
Fig. [B| are the window functions for the various configurations of the experiments 
considered. 

On the largest scales corresponding to small I, new COBE jl8| and Tenerife 
|89t results improve the power spectrum normalisation, whilst significant gains 
in knowledge at high I are provided by new results from the Saskatoon and CAT 
experiments. The full data set spans a range of 2 to ~ 700 in I, sufficient to test 
for the main Doppler peak out to Q = 0.1. We take the reported CMB detections 
and convert them to a common framework of flat bandpower results Ji^, |24|] as 
given in Table |j| This is carried out as follows. 

The CMB anisotropy measurements are converted to bandpower estimates 
AT; ± a assuming in each case a flat spectrum of C; centred on the effective 
multipole l e (see below) of the window function, k and l u represent the lower 
and upper points at which the window of each configuration reaches half of its 
peak value. In order to use the observed anisotropy levels to place constraints 
on the CMB power spectrum one must in general know the form of the Ci under 
test. However, in most cases the form of Ci can be represented by a flat spectrum 
Ci oc C'2/(l(l + 1)) over the width of a given experimental window, so that the 
bandpower is ATJ/T = y/ C o bs(0)/I(Wi), where we define I(Wi) according to 

Bond IHIH as I( ~ W ') = IX 2 ( Z + a5 ) W VW / + 1 )- This bandpower estimate is 
centred on the effective multipole l e = I (IWi) / 1 (Wi) . In many instances exper- 
imenters now report results directly for a flat spectrum and when this is not so 
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Table 3. Details of data results used 
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Figure 6. The data points from Table || are shown compared to the best 
fit analytical CDM model. The dotted and dashed lines show the best fit 
models which are obtained when the Saskatoon calibration is adjusted by 
±14%. The data points from the MAX experiment are shown offset in I 
for clarity 



we have converted the quoted power in fluctuations into the equivalent flat band 
estimate. Each group has obtained limits on the intrinsic anisotropy level using a 
likelihood analysis (see e.g. Hancock et al. [|9l| ), which incorporates uncertai nties 
due to random erro rs, sampling variance [1169] and cosmic variance 166 
Section H and 



11671 (see 

158]). The errors in AT; quoted in column 3 of Table |3| are at 68 



% confidence and have been obtained by averaging the difference in the reported 
68% upper and lower limits and the best fit ATj . Since the form of the likelihood 
function is in general only an approximation to a Gaussian di strib ution this aver- 
aging introduces a small bias into the results (see Section ^j, [154 1). Results from 
the MSAM experiment are not included here, because they do not provide an 
independent measure of the power spectrum since their angular sensitivity and 
sky coverage are already incorporated within the Saskatoon measurements. (If 
the correlation is small, as recent information seems to indicate, then this ex- 
periment may be included in the analysis; this is done in Section ^|) Netterfield 
et al. J136] ] report good agreement between the MSAM double difference results 
and Saskatoon measurements, although the discrepancy with the MSAM single 
difference data is yet to be resolved. The data points from Table ^| are plotted 
in Figure |^, in which the horizontal bars represent the range of I contributing 
to each data point. There is a noticeable rise in the observed power spectrum 
at I ~ 200, followed by a fall at higher I, tracing out a clearly defined peak in 
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the spectrum. In the past several groups [ 17d| , 106, 147, frij have attempted to 
determine the presence of a Doppler peak, but only now are the data sufficient 
to make a first detection and to put constraints on the closure parameter SI. As 
a first step, we adopt a simple three parameter model of the power spectrum, 
which we find adequately accounts for the properties of the principal Doppler 
peak for both standard Cold Dark Matter (CDM) models Q and open Uni- 



verse (Si < 1) models 106]. The funct ional form chosen is a modified version of 



that used in Scott, Silk & White |170| — we choose the following 

1(1 + 1)C, = 6C 2 (l + - Apea " ) /(l+ , f pe °* ) (.",?) 



i + y(0V ' V i + y(2) 2 

where y(l) = (log 10 I — /o(/io(220/\/r2))/0.266. In this representation C2 speci- 
fies the power spectrum normalisation, whilst the first Doppler peak has height 
Ap ea k above C2, width log 10 / = 0.266 and for SI = 1.0 is centred at I ~ 220. By 
appropriately specifying the parameters C2, Apeak and SI it is possible to repro- 
duce to a good approximation the Ci spectra corresponding to standard models 
of structure formation with different values of SI, Sib and Hq. Such a form will 
not reproduce the structure of the secondary Doppler peaks, but we have checked 
the model against the overall form of the SI = 1 mod els of Efstathiou and the 



open models reported in Kamionkowsky et al. |106| and find that this form 
adequately reflects the properties of the main peak. This satisfies our present 
considerations since the current CMB data are not yet up to the task of discrimi- 
nating the secondary peaks. Varying the three model parameters in equation (^7]) 
we form Ci spectra corresponding to a range of cosmological models, which are 
then used in equation ([si]) to obtain a simulated observation for the ith experi- 
ment, before converting to the bandpower equivalent result ATJC2, A pea k, Si](i). 
The chi-squared for this set of parameters is given by 

2 (r , v (AT^ji) - AT; [C*2 , Apeak , SI] (i)) 2 

j=l 1 

for the nd data points in Table ^ and the relative likelihood function is formed 
according to L(C2, A pea k, SI) oc exp(— x 2 (C-2, A pea k, Sl)/2). We vary the power 
spectrum normalisation C2 within the 95 % limits for the COBE 4-year data ]l8| 
and consider A pea k in the range to 30 and values of the density parameter up 
to SI — 5. The data included in the fit are those from Table ^| which with the 
exception of Saskatoon include uncertainties in the overall calibration. There is 
a ±14% calibration error in the Saskatoon data, but since t he S askatoon points 



are not independent this will apply equally to all five points 1 136 ] . The likelihood 
function is evaluated for three cases: (i) that the calibration is correct, (ii) the 
calibration is the lowest allowed value and (iii) the calibration is the maximum 
allowed value. In each case the likelihood function is marginalised over C2 be- 
fore calculating limits on the remaining two parameters according to Bayesian 
integration with a uniform prior. 



5.2. Results and D 



iscussion 



In Fig. the likelihood function obtained from fitting the model Ci spectra to 
the data of Table |^ is shown plotted as a function of the amplitude and position 
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(f2) of the Doppler peak. The highly peaked nature of the likelihood function 
in Fig. [?] is good evidence for the presence of a Doppler peak localised in both 
position (ft) and amplitude. In Fig. |^ we also show the 1-D conditional likelihood 
curve for Q, obtained by cutting through the surface shown in the left hand side 
of Fig. Q at the best-fit value of A pea k- The best fit value of Q. is 0.7 with 
an allowed 68% range of 0.30 < Q, < 1.73. In Figure ^ the best fit model, 
represented by the solid line, is shown compared to the data points, assuming 
no error in the calibration of the Saskatoon observations. The chi-squared per 
degree of freedom for this model is 0.9, implying a good fit to the data. The 
peak lies at I — 263^<!,4 9 corresponding to a density parameter fl = 0.70^1; the 
height of the peak is A pea k = Q.Olj'g- The dashed and dotted lines show the best 
fit models (fi = 0.701°;^, A peak = 11.0 j£g and fi = 0.6S±^|, A peak = 6.5±f;g 
respectively) assuming that the Saskatoon observations lie at the upper and lower 
end of the permitted range in calibration error. These likelihood results using 
the analytic form for the C; and the results from a chi-squared goodness of fit 
analysis using exact models (see below) imply that independent of calibration 
uncertainties in the data, current CMB data are inconsistent with cosmological 
models with £1^0.3. 

The analytic approximation to the true Ci such as we use here, is a useful 
general tool, but as a detailed check we have also applied the chi-squared good- 
ness of fit test to actual COBE normalised C; models. A detailed comparison 
of the CMB data with these models is considered in Section [| For example, 
considering a range of CDM models with varying Q, in order to find the low- 
est Q compatible with the observations, we have considere d ex act models with 
H = 50kms _1 Mpc -1 and Vi b = 0.03 for = 0.1 - 0.5 |o|]. We find that 
Q = 0.5 is allowed, fl — 0.3 and below are completely ruled out (95% confidence) 
and fl — 0.4 is excluded unless all the Saskatoon points have the minimum al- 
lowed calibration. This result is confirmed in Section ^|for a more complet set of 
open models. 



5.3. Conclusions 

Our current results provide good evidence for the Doppler peak, verifying a crucial 
prediction of cosmological models and providing an interesting new measurement 
of fundamental cosmological parameters. This first estimate of the angular po- 
sition of the Doppler peak is used to place a new direct limit on the curvature 
of the Universe, corresponding to a density of SI = 0.7jlJ'4, consistent with a 
flat Universe. Very low density 'open' Universe models are inconsistent with this 
limit unless there is a significant contribution from a cosmological constant. In 
section ^ a detailed comparison of the CMB data is made with the theoretical 
power spectra predicted by a range of flat, tilted, open models and models with 
non-zero cosmological constant. 
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Figure 7. Top Figure: The likelihood surface for Q and A pea k- (The nomi- 
nal Saskatoon calibration is assumed. Bottom Figure: The 1-D conditional 
likelihood curve for 0. 
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6. Constraints on the Cosmological Parameters using current 
CMB observations-II 



In this Section we give a detailed comparison of the CMB data with the theoret- 
ical power spectra predicted by a range of flat, tilted, open models and models 
with non-zero cosmological constant. In this analysis we only consider scalar 
fluctuations and we do not include reionization. This is an extension of the work 
carried out in Section |^ in which we used an analytic approximation of current 
CMB models. Here we consider some of the exact fourth-year COBE normalized 
angular power spectrum with fourth-year COBE normalization. The precise form 
of the Doppler peaks depend on the nature of the dark matter, and the values of 
Q,, Qb, and Ho . Thus in order to use the medium-scale anisotropy results e.g. to 
provide additional leverage on the spectral index determination, it is necessary to 
adopt a given cosmological model. Each Section gives the results for each set of 
models of structure formation used in this intercom parison of data and models. 
For more information see also |89| p^ , 155, 157, 158]. Other related analysis usig 



a x 2 technique have been developped by, among others, J72j], [ 126 1 , etc. 

Ganga et al. (1997a, hereafter GRGS)) developed a technique to account for 
uncertainties, such as those in the beamwidth and the calibration, in likelihood 
analyses of CMB anisotropy data. This technique has been used in conjunction 
with theoretically-predicted CMB anisotropy spectra in analyses of the Gunder- 
sen et al. (1995) UCSB South Pole 1994 data (SP94), the Church et al. (1997) 
SuZIE data, the MAX 4+5 data (Tanaka et al. 1996; Lim et al. 1996), the Tucker 
et al. (1993) White Dish data, and the de Bernardis et al. (1994) ARGO data 
(GRGS; Ganga et al. 1997b, 1998; Ratra et al. 1998, 1999a, hereafter R99a). 
A combined analysis of all these data sets is presented in Ratra et al. (1999b, 
hereafter R99b) in which a joint likelihood analysis of these data sets is carried 
out. Bond & Jaffe (1997) have also analysed the SP94 data and the Saskatoon 
(Netterfield et al. 1997) data. A similar analysis of CMB anisotropy data from 
the Python I, II, and III observations performed at the South Pole (Dragovan et 
al. 1994, hereafter D94; Ruhl et al. 1995b, hereafter R95; Piatt et al. 1997) is 
described in |159| ] . The last Section gives a brief account of these other results in 
conjunction with the ones obtained here. 

6.1. Data 

We apply a fit procedure to three sets of CMB data first: the CMB data given in 
Table g of Section [|, where the COBE data consists in one multipole for COBE 
(namely C2, obtained from the value of the Q r ms- P s =18 [iK for the 4-year 
COBE data ^|) which we refer to as case (a); the same as the data of Table [5] 
apart from COBE where we consider 4 of the 8 multipo le b ands of the angular 
power spectrum extracted from the 4-year COBE data 



188], which we refer to 



as case (b), and the data consisting basically in the extension of data set (a) to 
include (1) the CAT1-97 point at I = 422 Q; (2) the inclusion of new points 
from experiments Python (Python III, Piatt et al 1997), MS AM (the 2nd and 
3rd flights, Cheng et al 1996, 1997, Ratra et al 1997), ARGO (Aries+Taurus 
region, Masi et al 1996), FIRS Q and BAM p90| , and with the latest cali- 



bration correction to the Saskatoon data (i.e. increased by 5%, Leitch, private 

41 



Table 4. Details of data results (c) used 
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communication); which we refer to as case (c) and is displayed in Table ^ and 
shown in Fig|^. In the case of data set (b), in order to avoid eventual correlations 
between these data points, we have applied the fitting analysis to the 4 even or 4 
odd multipole bands in conjuntion with the remaining data points. The results 
are similar whether we use the even or the odd multipoles. The use of 4 instead 
of 1 multipole is particularly useful in discriminating the value of the spectral 
index n. 

As already mentioned in Section ^| the different experi men t s sam ple different 
angular scales according to their window functions, Wi, §o|, p05j . The 4 win- 
dow functions corresponding to the COBE experiment were obtained considering 
Gaussian functions centere d in (I) with dispersion (A£) 2 , where the parameters 
are given in Tegmark 1996 [185]. 

The conversion of data to a common framework in terms of flat bandpower 
estimates [^, Q follows the same procedure as in Section ^ In Table [| we 
display these bandpower estimates AT; ± a where as before li and l u represent 
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Figure 8. The updated data set (c). 
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the lower and upper points at which the window function for each configuration 
reaches half of its peak value. As in Section [B|, the errors in AT; quoted in column 
3 are 68% confidence limits and have been obtained by averaging the difference 
in the reported 68% upper and lower limits and the best fit AT;. 



6.2. Method 

We pursued the same method as described in Section ^. The process by which 
the error bars quoted in Table ^ are computed introduces a small bias into the 
results since the likelihood function is in general only an approximation to a 
Gaussian distribution. In order to assess the significance of this we applied the 
fitting analysis to an alternative set of data which consisted of the mean value 
of the data and corresponding error bars. We find no significant difference in 
the results obtained. We consider that a given model offers an acceptable chi- 
squared fit when P{x 2 ) > 0.05. We increase the resolution of the parameters grid 
via interpolation procedures and compute the \ 2 and corresponding P(x 2 ) f° r 
each model. Including the 95% COBE normalization uncertainty and increasing 
the resolution of the parameter grid using interpolation procedures, we computed 
the best-fit model, its chi-squared value, the corresponding probability and the 
likelihood function. For each parameter both the conditional and marginal dis- 
tributions were obtained as corresponding confidence intervals. 



6.3. Scale invariant flat CDM models 

This work was initially applied to scale invariant (n — 1), flat (f2 = 1) CDM 
models with n b =0. 01-0.3 and # =30-100 km s" 1 Mpc" 1 (kindly provided by G. 
Efstathiou). We considered these models normalised to the 4- year COBE value 
of Qrms-ps =18 /iK and applied the analysis to the data tabulated in Table ^| 
of Section |B| (case (a)). In Fig. ^ we have considered these COBE normalised 
CDM models. A dot is placed in the appropriate place in the parameter space 
if the exact power spectrum corresponding to these parameters gives a fit to 
the data in Table ^| of Section |E] with an acceptable x 2 value (-P(x 2 ) > 0.05). 
A blank is left at that position if not. Overlying these power spectrum con- 
straints is the limit 0.009 < Q. b h 2 < 0.02 provided by nucleosynthesis of the 
light elements As shown, the models offering an acceptable chi-squared fit 
(-P(X 2 ) > 0.05) to the CMB power spectrum whilst simultaneously satisfying nu- 
cleosynthesis constraints Ji^ , are those with 0.05 < Q, b < 0.2, 30kms _1 Mpc -1 < 
Ho < 50kms _1 Mpc -1 . Considering the highest Saskatoon data calibration the 
constraints become 0.1 < fif, < 0.2, 30 km s" 1 Mpc" 1 < H < 35kms _1 Mpc" 1 . 
Allowing for the lowest Saskatoon data calibration relaxes the constraints up to 
H = 70 km s" 1 Mpc" 1 and 0.02 < Q b < 0.2 (see Hancock et al. j89|). In gen- 
eral, rece nt optica l and Sunyaev-Zel'dovich observations of the Hubble constant 



144 , |68| , |L07| , p.141 imply H in the range 50 — 80km s 1 Mpc 1 . Since the paper 



Hancock et a l. |8Sfl was first submitted, an alternative comparison of CMB data 
with models [ p.2q ] has appeared, which supports our conclusions that low values 
of Ho are favoured by the current CMB data. Alternatively we generated the an- 
gular power spectrum, Cj , for these models using a CMB code kindly provided 
by U. Seljak and M. Zaldarriaga. The grid considered has Hubble constant, Ho, 
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Figure 9. The COBE normalised CDM models with acceptable chi-squarcd 
fits to the CMB data (for case (i) calibration) are plotted as dots in the 
ilb — H space. Overlayed are the constraints on Qi,h 2 imposed by nucle- 
osynthesis. 



Ho Qrms — ps 

nBBN 

Best Model 30 0.02 17.64 

Conditional la) 30-50 0.02-0.15 16.02-19.44 
Marginal(lcr) 30-60 



BBN 

Best Model 30 0.1 16.38 

Conditional la) 30-40 0.1-0.19 14.94-17.82 

Marginal(lcr) 30-50 0.04-0.15 16.20-19.98 



Table 5. Results of the fitting analysis considering: one multipole for COBE 
experiment (a), that Saskatoon calibration is correct (i) and with (BBN), 
without (nBBN), BBN constraints for scale invariant flat CDM models. 



45 



model n (68% conf. interval,) 







(a),(i) 


(b),(i) 


(a),(i,ii,iii) 


(b),(i,ii,iii) 


Q b 


= 0.3 
= 0.2 


0.94 ±0.06 


0.99 ±0.07 


0.94 ±0.1 


99+ 11 


h 
Q b 


= 0.45 
= 0.1 


1.03 ±0.07 


1.08 ±0.07 


1.03 ±0.1 


1.08 ±0.1 


ho 


= 0.5 
= 0.07 


1.08 ±0.07 


1.11 ±0.07 


i.08l° :S 


i -I I +0.09 
i - L± -0.11 



Table 6. Results of the fitting analysis for tilted flat models for: (a) data 
with 1 COBE multipole, (b) data with 4 COBE multipoles; (i) Saskatoon 
calibration is correct and (ii) lowest (iii) higher Saskatoon calibration. 



values in the range Hq = 30 — 80 kms" 1 Mpc _1 in steps of 5 kms -1 Mpc -1 , baryon 
density, £l b , values in the range £l b = 0.02 — 0.22 in steps of 0.01 and f2o = 1. We 
applied the same analysis to these models considering the data set (a), assuming 
that the Saskatoon calibration is correct (i). The results are displayed in Table H. 
In order to set constraints on the parameters of these set of models we consider 
the 68% c.i. for the marginalized distribution of each parameter constrained to 
the BBN limits, these encompass 30 kms -1 Mpc -1 < H < 50 kms -1 Mpc -1 , 
0.04 < fl b < 0.15 and 16.2 ^K< Q rms -ps < 19.98 /xK. It should be stressed 
that these results were obtained considering only the case (i) for the Saskatoon 
calibration. This analysis confirms the previous results obtained using a cruder 
method. 



64. Tilted flat CDM models 

Initially we considered some of the best candidates for flat CDM models: (1) 
h = 0.3, fl b = 0.2, (2) h = 0.45, Q b = 0.1 and (3) ho = 0.5, Q b = 0.07 and 
allowed for the tilting of the power spectrum to vary within n — 0.7 — 1.4. The 
fitting for the spectral index n of the primordial fluctuations was done considering 
the two sets of data (a) and (b) mentioned above while at the same time allowing 
for the 2cr COBE normalization uncertainty. The results are displayed in Table M. 
The third and fourth column display the uncertainties obtained considering the 
three calibration cases and taking extreme limits. Comparing e.g. first and second 
column we conclude that in general the effect of using 4 multipoles for COBE data 
produces an increase of n. In particular fixing Ho = 50 kms" 1 Mpc -1 we find 
n = 1.1±0.1 (68 % confidence interval) allowing for the three calibration cases and 
taking extreme limits. This tight limit rules out a significant contribution fr om a 



gravity wave background for these models, in the case of power law inflation [ 12C ] , 
but is consistent with the prediction of n ~ 1.0 for scalar fluctuations generated 
by inflation. In Fig. [l(] on the left hand side plot we present the model (3) and its 
68 % confidence interval assuming the case (i) for the calibration. On the right 
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Figure 10. Tilted model with h = 0.5, Sl b = 0.07, n = 1.1 and its 68 
%confidence interval (Left hand side Figure), likelihood function for the 
parameter n (Right hand side Figure). 



hand side plot is displayed the likelihood function for the parameter n for model 
(3) which peaks at n = 1.1 showing a significant estimate of the spectral index 
parameter based on the act ual CMB data. For details about the normalization 
of these models see e.g [158 . 

We then considered the 4-year COBE normalized tilted flat CDM models. 
We generated the angular power spectrum, C^ S \ using a CMB code provided by 
U. Seljak and M. Zaldarriaga. The grid considered has Ho values in the range 
Ho = 30 — 80 kms -1 Mpc -1 in steps of 5 kms" 1 Mpc~\ Sit values in the range 
Sib = 0.02 — 0.22 in steps of 0.01 and tilt in the range n = 0.6 — 1.4 in steps of 
0.1 for SI = 1. The results for case (a) are displayed in Table with and without 
superimposed BBN constraints as well as for the updated data set (c). 

Considering the BBN constraints we obtain for the marginalized distribution 
of the tilt a 68% c.i. of 0.85 < n < 1.18 while not including this constraint 
implies an n in the range 0.78-1.12 (68 %). Adding new data points to data set 
(a) to get data set (c) affects the results by decreasing the value of the tilt n from 
0.92 to 0.88. The marginal distribution suffers a slight variation with an interval 
0.85-1.18 for case (a) and 0.80-1.15 for case (c). Both results are consistent with 
the prediction of n ~ 1.0 for scalar fluctuations generated by inflation. 

Fixing Ho = SOkms" 1 Mpc -1 the results obtained are tabulated in Table H. 
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Table 7. Results of the fitting analysis considering: one multipole for COBE 
experiment (a) and the updated data set (c), that Saskatoon calibration is 
correct (i) and with (BBN), without (nBBN), BBN constraints for tilted 
flat CDM models. 



The 68% c.i. for the marginalized distribution of the tilt suffers a shift towards 
higher values of n when we include the BBN constraints from 0.84-1.16 to 0.93- 
1.21. This is expected to be so because the BBN limits will not allow the higher 
values of Q.b which seem to fit better the data. 

In order to constrain Ho for tilted CDM models consider Table |^. The inclu- 
sion of BBN constraints reduces the 68% c.i. for the marginalized distribution 
from 30 km s _1 Mpc -1 - 70 km s" 1 Mpc" 1 to 30 kms" 1 Mpc" 1 - 55 km s" 1 Mpc" 1 
assuming that Saskatoon calibration is correct. Fixing n = 1 reduces the upper 
limit of the Ho c.i. in both situations. The marginalized distribution of the 
Hubble constant when considered in conjunction with BBN limits gives Ho in 
the range 30 kms" 1 Mpc" 1 - 55 kms" 1 Mpc" 1 for case (i) calibration. 

In Table |^ we also display the results obtained for 4 multipole bands for 
COBE data (b), assuming that Saskatoon calibration is correct and including 
the BBN constraints. The best fit value of n is n = l-15±g;i| (68 % c.i.) while 
the marginalized distribution gives n in the range 0.93-1.24 with 68 % confidence. 
The information contained in these 4 multipole bands seems to point for the ex- 
istence of more power in the smaller scales probed by COBE and for the decrease 
of the normalization at low I. The marginal distribution suffers an insignificant 
variation with an interval 0.93-1.21 for case (a) and 0.93-1.24 for case (b). These 
results are consistent with the inflationary prediction of n ~ 1. The presence of a 
gravity wave component in our models would require even larger values of n than 
those derived above. Considering Section [l| this reaffirms our conclusion that a 
significant gravity wave contribution is unlikely. 
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Qrms—ps 


(a),nBBN 
Best Model 
Conditional(l(j) 
Marginal(l(j) 


0.22 
0.15-0.22 


0.91 
0.81-1.0 
0.84-1.16 


17.68 
15.89-19.66 


(a),BBN 
Best Model 
Conditional(lcr) 
Marginal(l(j) 


0.08 
0.05-0.08 


1.07 
0.98-1.18 
0.93-1.21 


16.79 
15.14-18.82 


(b),BBN 
Best Model 
Conditional(lcr) 
Marginal(lcr) 


0.08 
0.05-0.08 


1.15 
1.03-1.26 
0.93-1.24 


14.20 
12.67-15.73 



Table 8. Results of the fitting analysis for Ho = 50kms 1 Mpc 1 for case 
(a) and (b) case (i) calibration with (BBN) and without (nBBN) BBN 
constraints for tilted flat CDM models. 



6.5. Open models 

We initially had access to the open models (fi < 1) with values for (Q,ho) (where 
ho = I/oAlOOkms^Mpc- 1 )) of (0.1,0.75), (0.2,0.7), (0.3,0.65), (0.4,0.65), 
(0.5,0.6) kindly provided by N. Sugiyama. Applying the analysis for case (b) 
ie 4 multipole bands for the COBE experiment, we conclude that Q < 0.3 is 
not compatible with the data, while the value of Q = 0.4 is only allowed for 
the lower Saskatoon calibration data. A value of fi = 0.5 offers an acceptable 
chi-squared fitting for both case (i) and (ii) calibrations. This confirms the re- 
sults obtained using the approximate formula of Section |E] reassuring us that 
the generalized parametrized formula previously used constitutes a good approx- 
imation given the uncertainties in the data. We have also considered a more 
com plet e set of open models with A = and open-bubble inflation spectrum 
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106 : the grid considered has h values of 0.3, 0.5, 0.6, 0.7 and 0.8 (where 
h~= 717/(100 kms" 1 Mpc" 1 )), and baryon density Q b values of 0.01, 0.03, 0.06, 
together with Qb = 0.0125/i -2 and 0.024/i -2 for each of the above values of h. 
The Q range considered is 0.1 to 1.0 in steps of 0.1. (These models were kindly 
provided by N. Sugiyama). The unavaiability of exact models for fi > 1 limits 
some of the statistical conclusions we can draw here, but the results are still of 
interest. Assuming case (i) for the calibration, data set (a), and allowing the 
model normalisation to vary within the two sigma COBE limit we find that the 
best fit model has Q = 0.7, H = 50 kms -1 Mpc" 1 and Q b = 0.096. This best 
fit value of Q gives good agreement with the results obtained using the analytic 
approximation. Marginalizing over the other parameters, we obtain an allowed 
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Table 9. Results of the fitting analysis for data set (a) and (c), and case 
(i) calibration for open and flat models. 



68% range for Q of 0.5 < Q < 1.0. (The upper limit of 1 is due to the cutoff in 
the range of models considered.) The results are presented in table assuming 
a similar amplitude for the the 95% uncertainty normalization as for the flat 
models. The confidence intervals for the marginal distributions of the parame- 
ters seem to be insensitive to the inclusion of new data points (as in (c)) and in 
particular of data from MSAM experiment. However, for data set (c), the best 
fit value of Hq increases and of Qb decreases slightly, when compared with the 
results for case (a). For the updated data set (c) the best fit model has Q — 0.7, 
Ho = 60kms _1 Mpc -1 and Qb = 0.067. Fig. [ill shows the best fit open model 
for the fitting analysis applied to case (a) and (c) and assuming that Saskatoon 
calibration is correct (i) (left hand side plot); and the marginal distribution for 
the parameter £1 (right hand side plot). 



6. 6. Tilted flat A models 

We considered, in the first instance, some of the cosmological constant, A, mod- 
els with no reionization, using the Berkeley CMB code. The initially used 
model parameters were (A,h ,£l b ) of (0.6,0.5,0.05), (0.6,0.5,0.03), (0.6,0.8,0.06), 
(0.7,0.5,0.05), (0.2,0.5,0.05). These models seem to be consistent with the data, 
with large values of A providing the best fit. A models compared with a flat 
CDM model with ho = 0.45 and Qb = 0.1 assuming case (i) calibration. We 
then considered a more complete set of flat models with a contribution from a 
cosmological constant, A, on the top of which we add some tilt. These models 
angular power spectra were generated using the CMB code provided by U.Seljac 
and M.Zaldarriaga. The grid considered has cosmological constant A density 
Qa in the range 0.3-0.7 in steps of 0.1, Hubble constant Ho in the range 30-80 
kms -1 Mpc -1 in steps of 5 kms -1 Mpc -1 , baryon density Qb in the range 0.02- 
0.22 in steps of 0.01 and tilt n in the range 0.6-1.4 in steps of 0.1. For details 



about the normalization of these models see e.g [158 1 . The fitting analysis was 
applied to case (a) for the COBE experiment, and case (c), case (i) calibration 
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Figure 11. For data set (a) (top figure) and for data set (c) (bottom 
figure): Left hand side plot shows the best-fit open model; Right hand side 
plot shows the marginal distribution for the parameter ft. 
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Table 10. Results of the fitting analysis considering data set (a) and the 
updated data set (c), and case (i) calibration with BBN constraints for 
tilted fiat A models. (Q) Qrms-ps, (Qcobe) Q r ms- V s = 18 jjK. 



and case (2) for the BBN constraints. The results are diplayed in Table |l0j, with 
inclusion of the effect of tilt, with n varying from 0.6 to 1.4 in steps of 0.01. 

Only fixing n = 1 and the COBE normalization we reobtain that larger 
values of A provide the best fit. The opposite is true if we include the COBE 
normalization uncertainty and allow the tilt to vary. In the former case the 68% 
confidence interval for the marginal distribution places A in the range 0.3-0.6, Ho 
in the range 40-65 kms -1 Mpc -1 while in the last case we get A in the range 
0.4-0.7, H in the range 30-55 kms" 1 Mpc" 1 , tilt in the range 0.8-1.1. These 
confidence intervals remain unaffected whether we consider data set (a) or data 
set (c). However the inclusion of new data points produces a decrease of the 
best fit value of fit and Qrms-ps- The conditional distribution is not particularly 
discriminatory for the A parameter and most of the models seem to offer good fit 
to the data. One must stress that the models considered only contemplate values 
of Ho and Sib satisfying BBN constraints. Fig. |l2] shows the best fit for tilted A 
models (left hand side plot) and the marginal distribution of the A parameter 
(right hand side plot). 
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Figure 12. For data set (a) (top figure) and for data set (c) (bottom 
figure): Left hand side plot shows the best-fit model for tilted A model; 
results constrained to BBN limits. Right hand side plot shows the marginal 
distribution for the A parameter. 
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6. 7. Topological defect models 



The sit uation for models in which structure formation is initiated by cosmic 
strings 1 1 30 , 142 is now more complex, since some of the predictions for the 
power spectra for strings have recently changed. 



We initially considered a string model (J. Magueijo, priv. comm., [13C] ) 
and texture models with (r2,/i ,fi () ) = (l,0.5,0.05), (1,0.7,0.05) (N. Turok, priv 
comm.). Models in which structure formation is initiated by cosmic strings |130| 
are consistent with the data, but are excluded for the higher calibration Saskatoon 
data. Texture models are also only allowed if one considers the lower calibr ation 
for the Saskatoon data. Previous calculations for the cosmic strings model 1 130 ] 
and low Q CDM models both have the first Doppler peak occurring in roughly 
the same position in I, so it might be thought surprising that only the latter 
are eliminated by the current data. This is traceable to the form of the low 
Q power spectra in the range I ~ 10-100, where in order to match the COBE 
normalization at low I, the models are forced to have values which are too low 
compared to the data over the intermediate angular scale range. It is likely, 
however, that these string models will be eliminated if more accur ate d ata on 
the CAT range of angular scales confirms the existing CAT results [171|. More 
recent calculations of t opo logical defect theories indicate that the Doppler Peak 
is strongly suppressed [ 142 1 and these predictions are likely to be ruled out by the 
current data. The latest cosmic string power spectrum predictions now include 
a cosmological constant (Avelino et al 1997, 1998, Battye et al 1998). These 
predictions succeed in recovering a significant 'Doppler peak' in the power spec- 
trum, but have peak power for I in the range 500 to 600, at variance with the 
trend of the current experimental results. In Fig. jl3] we present the data points 
of data set (b) compared with a cosmic string model, the best candidate of the 
texture models (with f2 = 1, ho = 0.5 and Qb = 0.05) and a flat CDM model 
with ho = 0.45 and fij, = 0.1. We emphasize that the topological defect theories 
shown in this plot represent the initial calculations of the Ci's and do not include 
the most recent predictions for the power spectrum, for these models. On the 
top plot we present the same data points compared now with candidates of both 
the open and A models and with the same flat CDM model. As shown in Fig. [Hj] 
the peak of the amplitude of the angular power spectrum for these models shifts 
towards higher multipole I (i.e. towards smaller angular scales) with respect to 
the standard CDM model. The position of the first Doppler peak for A models 
is ~ that of the standard CDM model as mentioned in Section [| 



6.8. Conclusions 

We have presented an intercomparison of the CMB data and the theoretical 
fourth-year COBE normalised power spectra predicted by several models. For 
flat, scale invariant, Cold Dark Matter dominated Universe, we used our re- 
sults in conjunction with Big Bang nucleosynthesis constraints to determine the 
value of the Hubble constant as Ho = 30 — 70kms~ 1 Mpc -1 for baryon fraction 
fib = 0.02 to 0.2. We also considered a more detailed analysis of the fitting 
via calculation of the probability distribution function of the parameters. We 
obtained for the marginalized distribution of the parameters constrained to the 
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Figure 13. The data points (b) compared to the exact forms of the Ci 
for on the top figure: an fl = 1, ho = 0.45 and fib = 0.1 standard CDM 
model (bold line), an fi = 0.5, h = 0.6 open model (dot line) and a 
flat f^A = 0.7, ho = 0.5 and fif, = 0.05 model (dashed line). For on the 
bottom figure: the same standard CDM model compared with a cosmic 
string model (dashed line) and an = 1, ho = 0.5 and = 0.05 texture 
model (dot line) . We emphasize that the topological defect theories shown 
in this plot represent the initial calculations of the Ci 's and do not include 
the most recent predictions for the f?<§wer spectrum, for these models. 



BBN limits, considering that Saskatoon calibration is correct, a 68% c.i. of 
30kms -1 Mpc -1 < H < 50kms -1 Mpc -1 , 0.04 < fl b < 0.15 and 16.2 fiK< 
Qrms-ps < 19.98 fiK, confirming previous results. Tilting the initial fluctuations 
away from scale invariant and fitting to the AT; we can delimit, n, using the 
two data sets: one COBE multipole (a) and four COBE multipole bands (b) for 
the same range in I. Assuming e.g. a CDM model without gravity waves with 
Ho = 50kms -1 Mpc -1 and fij, = 0.07 we find the primordial spectral index of 
the fluctuations to be n — 1.1±0.1 (68% c.i.), with best fit values of Qrms-ps =15 
fj,K for n = 1.1 (allowing for the three calibration cases). This COBE normalised 
tilted model predicts a bias of the order b ~ 0.7. So although this model fits 
well the CMB data alone it may not give a realistic scenario when we consider 
jointly the CMB data and the observed galaxy clustering. Aplying a more precise 
method to models with Ho = 50 km s -1 Mpc -1 , we found for the marginalized 
distribution of n a 68% c.i. of 0.93 < n < 1.21 for case (a) and 0.93 < n < 1.24 
for case (b), assuming that Saskatoon calibration is correct. Including all our 
set of Ho values we get for the marginalized distribution of the tilt a 68% c.i. 
0.85 < n < 1.18 for case (a) and 0.80 < n < 1.15 for the updated data set (c). 
These results are consistent with those from large scales alone given in Section kL 
and are in close agreement with the inflationary prediction of n ~ 1.0. They also 
state that a significant gravity wave contribution is unlikely. 

Including open models in our analysis we obtain results consistent with 
those of Section ^| For data set (a) the best fit model has Q — 0.7, Ho = 
50kms -1 Mpc -1 , Qb = 0.096 and Qrms-ps =15.24 jj,K while for case (c) the 
best fit model has Q = 0.7, Ho = 60kms -1 Mpc -1 , fit = 0.076 and Q r ms- P s 
=15.45 /iK. Marginalizing over the other parameters we obtain a 68% c.i. for 
f2 of 0.5 < Q < 1.0, where the upper limit of 1 is due to the range of models 
considered. For tilted models with non zero cosmological constant we obtain that 
lower values of A provide the best fit. For data set (a) the best fit has A=0.3, 
Ho = 35kms -1 Mpc -1 , fit = 0.16, n — 0.9 and Qrms-ps =18 /iK, which changes 
slightly for case (c) for the values of fit — 0.12 and Qrms-ps =17.82 /iK. The 
marginal distribution places A in the range 0.4 < A< 0.7 (68%) and the tilt in 
the range 0.8 < n < 1.1. In common with other recent work on fitting to current 
CMB data (e.g. Lineweaver & Barbosa 1998) a tendency to low Hq values (except 
in the case of open models) is found, although it is clear from the marginalised 
ranges that the statistical significance of this is not yet high. Work has been 
developed to constrain cosmological parameters using combined data from CMB 
and large scale structure observations. As an example, among others, see e.g. 
20C], etc. For comparison with the results of Webster el al (1998), who worked 
jointly with recent CMB and IRAS large-scale structure results, we note that a 
flat A model with normalization fixed to the COBE results and n = 1 yielded 
a best fit of f?A = 0.7 and h = 0.6. In conjunction with the CAT points, new 
results from the Python and QMAP experiments may soon be very significant in 
ruling out some models and in further delimiting the Doppler peak, sharpening 
up these parameter estimates. Recent OVRO results at t ~ 590 (Leitch, private 
communication) agree well with the values found by CAT. The joint CAT and 
OVRO results will clearly be very significant in constraining the latest cosmic 
string power spectrum predictions, which now include a cosmological constant 
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For comparison we give a summary of some of the results obtained from data 
analysis of other experiments and in particular the results from a joint likelihood 
analysis of combined data sets for Open and Flat-A CDM Cosmogonies. The full 
SP94 data set is most consistent with Qo ~ 0.1 — 0.2 open models and less so 
with old (to > 15 — 16Gyr), high baryon density (f2& > 0.0175fo -2 ), low density 
(Oo ~ 0.2 — 0.4), fiat- A models. The SP94 do not rule out any of the models 
considered at 2a level j74j. The SP94 data subsets model normalizations are 
mostly consistent with those deduced from the 2 yr COBE-DMR data. The 
same analysis now applied to MAX data give results consistent with those drawn 
from the SP94, mentioned above, but only for open models with fio ~ 0.4—0.5 the 
MAX and DMR normalizations overlap, the same happens for f2o = 1 and some 
higher h, lower fij,, low-density flat-A models [j75[. The White Dish data does not 
give 2a detections of anisotropy for the models tested. The 2a upper limit for 
the flat bandpower angular power spectrum is 150 ^K(<3 rms _ ps < 96 fiK) (HPD 
prescription, see |7^] for further details). The upper limits for all models tested 
are consistent with those derived from the COBE-DMR data [149|. The same 
happens for the SuZIE data with a 2a upper limit for the fiat bandpower angular 
power spectrum of 40 /iK(Q rms _ ps < 26 /iK) The ARGO data favors 

an open (spatially-flat- A) model with Qq = 0.23(0.1). At the 2a confidence 
level model normalizations deduced from the ARGO data are consistent with 
those drawn from the SP94, MAX 4+5, White Dish, and SuZIE data sets. The 
Argo open model normalizations are also consistent with those deduced from the 
DMR data. However, for most spatially-fl at-A models the DMR normalizations 
are more than 2a above the ARGO ones 1 148 ] . The joint likelihood analysis of 
SP94, ARGO, MAX, White Dish, and SuZIE data favors an Q ~ 1(1) open 
(flat-A) model. Excluding SuZIE data, an Qo — 0.3(1) open (flat-A) model 
is favored. Considering only multi-frequency data with error bars consistent 
with sample variance and noise considerations an Oo — 0.1(1) open (flat- A) 
model is favored. However, the data do not rule out other values of Qq in the 
flat-A models and other values of Q,bh 2 in both models. The small-scale data 
normalizations are consistent at 2a level with those derived from the DMR data. 
The Python 1,11 and III data combined favors an open (spatially-flat-A) model 
with flo ~ 0.2 (0.1). At the 2 a confidence level model normalizations deduced 
from the combined Python data are mostly consistent with those drawn from the 
DMR, UCSB South Pole 1994, ARGO, MAX 4 and 5, White Dish, and SuZIE 
data sets. While the statistical significance of the constraints on cosmological 
parameters is not high, it is interesting that the combined Python data favor 
low-density, low Qb h 2 , young models, consistent with some of the indications 
from the combinations of CMB anisotropy data considered by R99b. The x 2 
analysis presented here favour high-density open models in variance with R99b, 
though in the latter case the results favour an Q, — 1 model when the smallest 
angular scale SuZIE data is included in the joint analysis. In our analysis here we 
include experiments such as Saskatoon and CAT which also probe higher values 
of I when compared with the remaining observations considered. The difference 
between these two analysis might be due to the fact that they use data from 
different combinations of experiments. It also seems to indicate that inclusion of 
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observations probing higher values of the multipole I, which do better pin down 
the position of the Doppler peak, shifts the density parameter to a larger value. 
However, given the large error bars, both analysis are consistent for f2 ~ 0.4 — 0.5 
at 2a level. 



7. Conclusions 

This paper has presented an intercomparison of CMB anisotropy experiments 
and models of structure formation. On large angular scales the CMB fluctua- 
tions may probe the initial power spectrum of the fluctuations, and in particular 
may be used to constrain the normalization and the spectral index of an initial 
power-law power spectrum. Results for both parameters were obtained for the 
Tenerife experiment, and an improved result achieved from the combined analy- 
sis of COBE and Tenerife data. These experiments in conjunction with medium 
angular scale experiments cover angular scales ranging from 10° to 0.25°, probing 
the shape of the CMB 'angular power spectrum' all way up to I ~ 700, offering an 
unique opportunity of constraining model parameters. We carried out with a full 
comparison of these recent detections of CMB anisotropies in order to constrain 
e.g. Oo, Ho, fit, n, A etc. A summary of the results from each of these areas is 
now given. 

An initial intercomparison of some CDM models with 'old' CMB data (avail- 
able at the time the work was carried out) ruled out isocurvature models. 

A summary of the analysis of the cosmic microwave background structure 
in the Tenerife Dec=+40° data has been presented. Our analysis demonstrates 
the existence of common structure in independent data scans at 15 and 33 GHz. 
For the case of fluctuations described by a Gaussian auto-correlation function, a 
likelihood analysis of our combined results at 15 and 33 GHz implies an intrin- 
sic rins fluctuation level of 48l^g /iK on a coherence scale of 4°; the equivalent 
analysis for a Harrison-Zel'dovich model gives a power spectrum normalisation of 
Qrms-ps = 22tg° A*K. The fluctuation amplitude is seen to be consistent at the 68 
% confidence level with that reported for the COBE two-year data for primordial 
fluctuations described by a power law model with a spectral index in the range 
1.0 < n < 1.6. This limit favours the large scale CMB anisotropy being domi- 
nated by scalar fluctuations rather than tensor modes from a gravitational wave 
background, if power law inflation is right. The large scale Tenerife and COBE re- 
sults are considered in conjunction with observational results from medium scale 
experiments in order to place improved limits on the fluctuation spectral index; 
we find n = 1.10 ± 0.10 assuming standard CDM with Ho = 50kms _1 Mpc -1 . 
For an improved grid of the parameters in conjunction with Big Bang nucleosyn- 
thesis constraints, we find 0.80 < n < 1.15 after marginalizing the distribution 
function over Ho , Qb, and Q rm s- V s ■ This constraint on the spectral index is in 
close agreement with the inflationary prediction of n ~ 1.0. 

A key prediction of cosmological theories for the origin and evolution of struc- 
ture in the Universe is the existence of a 'Doppler peak' in the angular power 
spectrum of cosmic microwave background (CMB) fluctuations. We presented 
new results from a study of recent CMB observations which provide the first 
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strong evidence for the existence of a 'Doppler Peak' localised in both angular 
scale and amplitude. This first estimate of the angular position of the peak is used 
to place a new direct limit on the curvature of the Universe, corresponding to a 
density of 57 = O.7IJ4, consistent with a flat Universe, making use of an analytic 
approximate expression to describe these models. This was further confirmed by 
the analysis conducted for the exact form of the angular power spectrum. Very 
low density 'open' Universe models are inconsistent with this limit unless there is 
a significant contribution from a cosmological constant. Marginalizing over the 
other parameters, we obtain an allowed la range for f2 of 0.5 < Q < 1.0. (The 
upper limit of 1 is due to the cutoff in the range of models considered.) The joint 
likelihood analysis carried out by R99b seems to favour low-density models in 
variance with the analysis carried out in Section ^. The difference between these 
two analysis might be due to the fact that they use data from different com- 
binations of experiments. In particular the \ 2 analysis presented here include 
experiments probing a region of high multipole I, which do better pin down the 
position of the Doppler peak. The effect of including SuZIE experiment in the 
analysis carried out by R99b is to shift the parameter density to a larger value, 
this reinforces our previous statement. However, given the large error bars, both 
analysis are consistent for Q ~ 0.4 — 0.5 at 2a level. For a flat Cold Dark Mat- 
ter dominated Universe with tilt we used our results in conjunction with Big 
Bang nucleosynthesis constraints to determine the value of the Hubble constant 
as Ho = 30 — 55kms _1 Mpc -1 (assuming that Saskatoon calibration is correct). 
The best fit A model has a small S1a content, the la confidence interval for the 
marginal distribution places A in the range 0.4-0.7. The latest cosmic string 
power spectrum predictions which now include a cosmological constant might be 
in variance with the set of the current experimental results. 

Our attempts to place increasingly more accurate limits on fundamental cos- 
mological parameters will undoubtedly place increasing demands on observers 
for ever improved accuracy until the point where the intrinsic cosmic variance 
becomes the dominant form of error. Confidence in the results to this tolerance 
level will probably require close co-operation between observers and the combi- 
nation of results from space-based, balloon-based and ground-based telescopes 
working over a range of frequencies. Potentially the most powerful observations 
will result from mapping overlapping/interlocking CMB fields with independent 
multi-frequency instruments. In contrast to the statistical results currently being 
reported, this method builds in the necessary redundancy to reduce systematic 
errors. As an example observations are currently in progress with the Tenerife 
instruments, which in conjunction with the COBE 4 year data should provide a 
high signal to noise map of the last scattering surface thus providing a two dimen- 
sional representation of the seed structures as compared to the simple 1-D scans 
reported here. Being on scales greater than the horizon size at recombination 
the form of such a map would reflect the structures generated in an inflationary 
driven phase in the very early universe. 

The existence of the Doppler peak has important consequences for the future 
of CMB astronomy, implying that our basic theory is correct and that improving 
our constraints on cosmological parameters is simply a matter of improved in- 
strumental sensitivity and ability to separate out foregrounds. New instruments 



59 



such as VSA [ 1 14 1 , MAP and the Planck Surveyor satellite [ 13 1 1 will provide 
this improved sensitivity and should delimit Q and other parameters with un- 
precedented precision. It is difficult using the CMB data alone to disentangle the 
variation due to simultaneous changes of different cosmological parameters, until 
such high resolution and high sensitivity maps of the sky are available. We can 
however use the information from other areas of astrophysics, for example the 
measurements of Large-Scale Structure (LSS). The CMB observations and the 
LSS measurements provide information about the power spectrum from comple- 
mentary regions, with some overlap, at different cosmological times. Since these 
two fields provide a larger coverage of the power spectrum, it is expected that 
their combined use will result in more severe constraints on the cosmological pa- 



rameters |17C[ . In order to include overlapping sky fields and correlated data 
points an improvement over the \ 2 analysis presented here is required. A % 2 
analysis is no longer suitable and this must be substituted by a correlation func- 
tion. However this may require the use of supercomputers and data compression 
techniques since matrix sizes in the likelihood analysis will be large. 
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